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APRESENTAÇÃO 


A coletânea “Física: O Universo e os Fenômenos Naturais” é uma obra que tem 
como foco principal a discussão científica por intermédio dos trabalhos que compõem 
seus capítulos. O volume abordará, de forma categorizada e interdisciplinar, resultados 
de pesquisas, relatos de casos e/ou revisões que transitam no pluralismo conceitual e 
epistemológico da Física e seu ensino. 

O objetivo central do livro é apresentar, de forma categorizada e clara, estudos 
desenvolvidos em diversas instituições de ensino e pesquisa do Brasil. A Física é uma 
ciência natural baseada em experimentos, medições e análises matemáticas com o 
propósito de encontrar leis físicas quantitativas para tudo, desde o nano mundo do 
microcosmo aos planetas, sistemas solares e galáxias que ocupam o macrocosmo. Temas 
diversos e interessantes são, deste modo, discutidos aqui com a proposta de fundamentar 
o conhecimento de acadêmicos, mestres, doutores e todos aqueles que de alguma forma 
se interessam pela Física e seus processos de ensino e de aprendizagem. 

Os autores do primeiro capítulo da obra, por meio dos parâmetros estruturais de um 
poço quântico de GaAs com barreiras de AlGaAs, simularam os níveis de energia, funções 
de onda e densidade de probabilidade, por meio de uma solução numérica da Equação de 
Schródinger, independente do tempo, para um poço finito. Com os resultados da interação 
elétron-buraco para os diferentes níveis energéticos do poço, o referido trabalho apresenta 
um exemplo direto e simples na solução de poços de potenciais quânticos reais. O segundo 
capítulo apresenta um trabalho experimental, no qual os autores utilizaram um termômetro 
de infravermelho de baixo custo para estimar a água precipitável na região de Florianópolis- 
SC. Fundamentado na dinâmica molecular clássica, os autores do terceiro capítulo discutem 
as propriedades termodinâmicas em sistemas binários e ternários compostos por óleo leve, 
C0O2 e salmoura para aplicações envolvendo recuperação aprimorada de petróleo. O quarto 
capítulo apresenta resultados de um estudo dos movimentos orbitais de detritos espaciais 
na vizinhança da Estação Espacial Internacional com o intuito de rastreá-los, visando evitar 
colisões entre eles e satélites artificiais operacionais. 

Os demais trabalhos apresentam instrumentos e metodologias para o ensino de 
Física. O quinto capítulo, por exemplo, usam a abordagem da modelagem científica de 
Bunge, adotando um objeto-modelo experimental e o incorporando na teoria da Mecânica 
Newtoniana para o desenvolvimento de um modelo teórico-prático validado empiricamente. 
Seguindo a perspectiva de propostas experimentais para o ensino de Física, o sexto capítulo 
discute a possibilidade de demonstração do efeito magnético de um imã em um solenoide 
por meio da construção de um Trem Magnético. O sétimo capítulo parte da grade comum 
dos conteúdos ensinados no curso de graduação em Física relacionados à dinâmica de 
rotações, para explicar o funcionamento do giroscópio, pião, tip-top e spinner, os quais, 
apesar de serem normalmente utilizados para recreação, possuem também aplicações 


educacionais e outras utilidades muito interessantes de serem exploradas na área de 
Ensino de Física. 

Além de práticas experimentais “manuais”, alguns trabalhos contemplam propostas 
para a ludicidade do ensino de Física. O oitavo trabalho incorporou o estudo de trajetórias 
bidimensionais de forma integrada, aliando a confecção de uma catapulta caseira em 
escala, no estilo Trebuchet, à realização de um estudo incisivo dos conceitos físicos 
necessários para descrição de seu funcionamento. A constituição de dados experimentais 
e simulações das trajetórias, foram obtidas por intermédio do software livre Tracker. O nono 
e último trabalho apresenta a potencialidade de recursos, como jogos educativos, como 
ferramenta de inclusão de alunos surdos, no que diz respeito aos processos de ensino e 
de aprendizagem de Física 1. 

Deste modo, essa leitura proporcionará um repertório de trabalhos bem 
fundamentados e com resultados práticos, obtidos por diversos professores e acadêmicos 
que arduamente desenvolveram seus trabalhos que aqui serão apresentados de maneira 
concisa e didática. Sabemos o quão importante é a divulgação científica, por isso 
evidenciamos também a estrutura da Atena Editora capaz de oferecer uma plataforma 
consolidada e confiável para estes pesquisadores exporem e divulguem seus resultados. 


Edson Ribeiro de Britto de Almeida Junior 
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RESUMO: No presente trabalho utilizamos 
parâmetros estruturais conhecidos de um poço 


Física: O Universo e os Fenômenos Naturais 


quântico de GaAs com barreiras de AlGaAs para 
simular os níveis de energia, funções de onda 
e densidade de probabilidade, através de uma 
solução numérica da Equação de Schródinger 
independente do tempo para um poço finito. Apartir 
dos resultados da simulação foi possível calcular 
a energia de transição elétron-buraco entre o 
nível fundamental e os primeiros níveis excitados 
do poço. Com os resultados da interação elétron- 
buraco para os diferentes níveis energéticos do 
poço, foi possível confrontá-los com um espectro 
de fotoluminescência da amostra obtido à 8 K, 
com o objetivo de identificar se o pico principal 
do espectro tem origem na transição do nível 
fundamental, primeiro nível excitado ou segundo 
nível excitado do poço. Portanto, este trabalho 
apresenta um exemplo direto e simples na 
solução de poços de potenciais quânticos reais. 
PALAVRAS  - CHAVE: Poço-quântico, 
Fotoluminescência, Simulação Numérica. 


OPTICAL ANALYSIS AND NUMERICAL 
SIMULATION OF THE ELECTRON-HOLE 
INTERATION IN À GAAS QUANTUM WELL 
WITH ALGAAS BARRIERS 


ABSTRACT: In the present work, we used known 
parameters of a GaAs quantum well with AlGaAs 
barriers to simulate energy levels, wave functions 
and probability density, using a numerical solution 
of the time-independent Schródinger equation 
for a finite quantum well. From the results of 
the simulation, it was possible to calculate the 
electron-hole transition energy between the 
fundamental level and the first excited levels of 
the well. With the results of the electron-hole 
interaction for the different energy levels of the 
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well, it was possible to compare them with sample's photoluminescence spectrum obtained at 
8 K, in order to identify whether the main peak of the photoluminescence spectrum is related 
to the transition from the fundamental level, the first excited level or the second excited level of 
the well. Therefore, this work presents a direct and simple example in the solution of potentials 
of real quantum wells. 

KEYWORDS: Quantum Well, Photoluminescence, Numerical Simulation. 


11 INTRODUÇÃO 


Nos dias atuais o estudo dos materiais semicondutores tem grande importância para 
o avanço da eletrônica e microeletrônica (DIAS, 1991). Dispositivos podem ser construídos 
através de diferentes estruturas e combinações de ligas semicondutoras, podendo formar 
estruturas binárias, ternárias e até mesmo quartanárias (DIAS, 1991). Uma dessas 
estruturas é a chamada poço quântico, onde se tem um elétron confinado entre barreiras 
de potencial finitas ou infinitas. A construção real desse sistema consiste em produzir um 
material semicondutor com um determinado gap de energia entre duas camadas de outro 
material com um gap maior (MORAIS et al., 2009). 

Conforme as características do poço quântico como largura e composição dos 
materiais que constituem as barreiras e o poço, os poços quânticos apresentam níveis 
quantizados de energia característicos para elétrons e buracos. Desta forma pode-se 
produzir diferentes sistemas com energias de transição desejadas, cuja aplicação maior 
reside na construção dos modernos e eficientes lasers de estado sólido amplamente 
empregados na indústria e meio científico. 

Utilizando parâmetros da amostra é possível realizar uma simulação teórica dos 
níveis de energia contidos na estrutura semicondutora, a fim de prever teoricamente o 
comportamento energético da mesma. Para realizar a simulação dos níveis de energia 
contidos em um poço quântico finito é necessário resolver a Equação de Schródinger, 
assim obtendo duas soluções que são utilizadas para obter de forma numérica os níveis 
de energia do poço. 

Neste trabalho apresentamos a revisão de um método simples para calcular 
teoricamente os níveis de energia em uma estrutura tipo poço quântico finito e real 
composto pelas ligas de Arseneto de Gálio (GaAs) com barreiras de Arseneto de Gálio 
Alumínio (Al Ga, As). Os valores obtidos teoricamente serão comparados com resultados 
experimentais obtidos através da técnica de espectroscopia de fotoluminescência. 


2 | MATERIAIS E MÉTODOS 


Bem tradicionais na mecânica quântica os poços de potenciais são soluções 
amplamente utilizadas, pois podem ter suas construções modeladas de várias formas para 
atingir objetivos tanto do ponto de vista acadêmico quanto em aplicações tecnológicas. O 
poço mais simples que podemos ter é um poço quadrado (simétrico em relação ao eixo 
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x) de potencial infinito. Uma aplicação simples do poço infinito é o confinamento de um 
elétron em uma placa metálica, muito fina (simulando uma situação 2D). Nesse sistema 
o confinamento acontece na direção da espessura da placa e o elétron pode se mover ao 
longo do plano (PAULINO et al., 2010). 

O poço de potencial infinito é um poço quântico onde sua principal característica 
é ter a energia potencial muito grande quando comparada a energia cinética da partícula 
(EISBERG; RESNICK, 1994), criando a possibilidade de confinar espacialmente um 
portador de carga em uma direção espacial (direção de crescimento), ao mesmo tempo em 
que o portador tem liberdade de movimento no plano perpendicular à direção onde ocorre 
o confinamento. Do ponto de vista da mecânica clássica, partículas confinadas neste poço 
de energia podem assumir qualquer energia, mas na mecânica quântica, apenas certos 
valores discretos de energia são possíveis (EISBERG; RESNICK, 1994). 

Na figura 1(a) podemos ver o esquema energético de um poço infinito. De —L/2 até 
L/2 temos um potencial igual a zero (região de dentro do poço com largura total L) e nas 
regiões fora do poço (de menos infinito até -L/2 e de L/2 até infinito) temos um potencial 
infinito (EISBERG; RESNICK, 1994). 


+00, x <-L/2 
VF(x)=4 0, -—L/2<xsL/2, (1) 
+es, x > L/2 
co Vo) co 
(a) (b) | 
+ V(X) 
Vo 
L/2 [0] L/2 -L/2 Õ L/2 


Figura 1: Representação esquemática (a) de um poço infinito e (b) de um poço finito. 
As energias e as funções de onda são encontradas exatamente ao resolver a 


equação de Schródinger independente do tempo, cujos resultados são dados pelas 
seguintes equações (EISBERG; RESNICK, 1994, ZETTILI, 2009): 
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n2hên? 


2mL? 


E, = n=1,23.. (2) 


E cos (=x); ni= 1,3,5 
Jêsm (=x);n=2,4,6.. 


Yn(X) = (3) 


onde À é a constante de Planck dividida por 2x, m é a massa do elétron, L a largura 
do poço e n é um número inteiro positivo que caracteriza o nível de energia. 

Nota-se que as funções de onda para n ímpar são simétricas, enquanto que para n 
par são anti-simétricas (ZETTILI, 2009). O módulo quadrado da função de onda fornece a 
densidade de probabilidade de encontrar uma partícula confinada em determinado estado 
(ZETTILI, 2009). 

Quando a energia potencial do poço não é suficientemente grande se comparada 
a energia cinética da partícula tem-se um poço quântico finito. O poço finito é um sistema 
quântico onde tem-se uma partícula confinada em um potencial (força atrativa exercida 
sobre a partícula) finito, em que o poço apresenta a característica de ser quadrado 
(EISBERG; RESNICK, 1994). Ao explorar esse problema bastante conhecido da mecânica 
quântica, pode-se observar e modelar fenômenos muito interessantes como, por exemplo, 
a quantização de energia do poço. Semicondutores podem ser utilizados para produzir 
estruturas reais com potencial do tipo poço quântico (FILHO; ALEXANDRINO; CABRAL, 
2018). 

Na figura 1(b) pode-se observar o esquema de um poço quântico quadrado finito. De 
menos infinito até — L/2 a energia potencial é V., de — L/2 até L/2 a energia potencial é zero 
e de L/2 até infinito o potencial é V, (EISBERG; RESNICK, 1994). 


Vo, x<-L/20ux> L/2 
V(x) = (4) 
0, —L/2 < x < L/2 


O ponto de partida para determinar os níveis de energia quantizados do problema de 
poço finito é a equação de Schródinger, para o nosso caso a mesma será unidimensional 
e independente do tempo (EISBERG; RESNICK, 1994). A equação de Schródinger é dada 


por: 


h2 d2p(x) 
zm dx? 





+ Voplo) = Ep(x). (9) 
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A solução geral para a equação (5), para a região dentro do poço (-L/2< x< L/2) é 
dada por (EISBERG; RESNICK, 1994): 


W(x) = Asen(k,x) + Bcos(k,x), (6) 


onde k, = /2mE/h. 


Para as regiões externas do poço, a solução da equação (5) é dada por (EISBERG; 
RESNICK, 1994): 


W(x) = Cekz* + DekzX, x<-L/2 (7) 
W(x) = FekzX + Gekzx | x>L/2 (8) 


onde kz = v2m(Vo — E)/h, e as letras de A até G são valores constantes, e k e k, 
são os correspondentes vetores de onda dentro do poço e das barreiras. 

A partir das equações (6), (7) e (8) e aplicando as condições de contorno 
(Vojx<-L/2oux> L/2 e 0,-L/2 <x < L/2) é possível obter duas soluções para o 
cálculo dos níveis de energia dentro de um poço finito (EISBERG; RESNICK, 1994). 














m'EL? mEL? m*(Vo— E)L2 
2h? Tan 2h? = 2h? (9) 
m*EL? m*EL? m'VoL?  m*EL? 
— | — = — — 1 
2hº Cot 2hº ) 2hº 2h? ( 0) 


A equação (9) descreve os níveis de energia pares (incluindo o nível fundamental), 
enquanto que a equação (10) descreve os níveis de energia ímpares (EISBERG; RESNICK, 
1994). 

Ambas as equações serão utilizadas posteriormente para efetuar o cálculo teórico 
e comparar os resultados com os obtidos experimentalmente em uma amostra de poço de 
GaaAs/AlGaAs. 


31 AMOSTRA 


A amostra analisada neste trabalho foi crescida através da técnica de MBE (Molecular 
Beam Epitaxy), sobre um substrato semi-isolante de GaAs (001) seguido por uma camada 
buffer de 1 ym de GaAs. Na sequência foi crescido uma super-rede com 30 repetições 
de [AI, (Ga, As/GaAs (50 À)], uma camada de 20 Á de AIAs e duas barreiras de 500 À 
de espessura de Al, Ga, , As confinando uma camada de GaAs de 150 À. Uma camada 


0,18 


0,18 
final de AIAs de 20 Á e uma Cap Layer de GaAs com 50 Á fecha a estrutura, conforme 


apresentado esquematicamente na figura 2. 


Física: O Universo e os Fenômenos Naturais Capítulo 1 Ei 


AlAs AlAs 


SL Alo 19 Gão g2ÃS Alo 19 Go g2ÃS 
É | | | 
A 
mm 


Figura 2 — Esquema estrutural da amostra analisada. 


Aamostrafoicaracterizadaatravés datécnicade espectroscopia de fotoluminescência, 
a temperatura de 8 K, com excitação proveniente de um laser de estado sólido emitindo em 
440 nme 6 uW de potência. A detecção foi feita utilizando um espectrômetro Andor de 19,3 
centímetros e uma CCD de Silício. 


41 RESULTADOS E DISCUSSÕES 

Para calcularmos os níveis de energia quantizados no poço vamos utilizar o modelo 
de poço quântico finito. Os gaps dos materiais que compõe o poço podem ser calculados 
em função da temperatura através das equações (BLAKEMORE, 1982): 


5,405 .1074 T? 


= oa 


(ev), (11) 


5,41.1074 T2 
Eg(T) = 1,519+ 1,155x + 0.37xº — “(T+204) | (eV), (12) 


A equação (11) é utilizada para calcular o gap do GaAs, enquanto para o cálculo 
do gap do Al Ga, As utilizamos a equação (12). Nestas equações T é a temperatura 
do material e o parâmetro x na equação (12) é a concentração de Al (alumínio) na liga 
ternária Al Ga, As. Utilizando a equação (11) com uma temperatura de 8 K a energia de 
gap do GaAs é 1,51884 eV e utilizando a equação (12), na mesma temperatura e com 
concentração de 18% de alumínio (x = 0,18) obtemos a energia de gap para o Al, Ga, AS 
igual a 1,/3872 ev. 

Neste tipo de estrutura o offset da banda de valência é determinado por E, = 0,17x 
(EXPUNOBI, 2002), onde x é a concentração de alumínio da estrutura (18%), assim 
temos E, = 0,0306eV. A partir dos valores de E, e dos gaps do GaAs e AlGaAs à 
temperatura de 8 K é possível determinar o offset da banda de condução (E,) como: 
Ec = 1,73872 — 1,51884 — 0,0306 = 0,18928eV. O offset da banda de condução equivale a 
profundidade do poço de potencial para os elétrons. De maneira similar é possível determinar 
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a profundidade do poço de potencial para os buracos na banda de valência através do 
offset da banda de valência (E,). A massa efetiva do elétron m* neste poço de GaAs é de 
0,067m, e para o buraco é de 0,53m,, onde m, é a massa do elétron (NAKWASKI,1995). A 
largura do poço para fins de cálculos será adotada como a largura nominal de crescimento 
da estrutura, 150 À. 

Os resultados dos cálculos teóricos para este poço quântico são apresentado na 


figura 3, com as indicações de Ev, Ec eos gaps do AlGaAs e GaAs. 


E>= 0,13796 eV 


E1= 0,06433 eV 


Ec = 0,18928eV 
gapaigass = 1,73872eV 


Eo = 0,016401 eV 


GaPGaas = 1,51884 eV 





Enn = 0,00217 eV 





Ev = 0,0306 eV 


150 À à 


Figura 3 — Representação esquemática dos níveis energéticos do poço quântico finito. 


Para obter teoricamente os valores dos níveis de energia do poço quântico (níveis 
de energia pares), foi utilizada a equação (9) resolvendo-a de forma numérica, onde m*é a 
massa efetiva do elétron (0,067m,), L a largura do poço quântico, h é a constante de Planck 
dividido por 27; Vo, a intensidade do potencial do poço quântico e E o nível de energia. 
Devido aos valores extremamente pequenos da constante de Planck, largura do poço e da 
intensidade do potencial é necessário reescrever estes parâmetros em termos de outras 
grandezas conhecidas, de modo a facilitar os cálculos numéricos. Assim a largura do poço 
foi reescrita em termos do raio do átomo de hidrogênio (ão = 0,53 À), a constante de Planck 
dividida por 27 foi normalizada (h =1), os offsets das bandas de valência e de condução que 
são equivalentes ao potencial do poço quântico foram convertidos em função da energia 
do átomo de Hidrogênio (£o = 27,2 eV, é duas vezes a energia do estado fundamental do 


átomo de hidrogênio) e a massa do elétron também é unitária (mo = 1), 
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Para a determinação dos níveis de energia na banda de condução e de valência 
foram utilizados os parâmetros da tabela 1. 


"Parâmetros Valores para banda de Valores para banda 
condução de valência 
L 283,01886 ao 283,01886 ao 
0,695882 .10"* £9 1,125.10"2 € 
m 0,067 mo 0,53 mo 


Tabela 1 — Parâmetros utilizados para determinação dos níveis de energia. 


Usando estes parâmetros e as equações (9) e (10) foi possível obter os resultados 
apresentados na tabela 2. 


CNN Ele) EM 
1º níveldopoço(Eo) —  0,000603 0,0416401 
2º nível do poço (Es) 0,002365 0,064330 
3º nível do poço (Ez) 0,005072 0,137960 
Nível da Banda de Valência (Ehn) 0,000079 0,002170 


Tabela 2 — Valores teóricos para os níveis de energia do poço (banda de condução) e para 
banda de valência em termos de E, e em unidades de ev. 


As energias E, E, e E, apresentadas na tabela 2 correspondem aos níveis 
quantizados na banda de condução, acessíveis para elétrons, enquanto que E,, é um nível 
quantizado da banda de valência, acessível para buracos presentes na região do poço 
quântico. 

A partir dos níveis de energia demonstrados na figura 3 juntamente com a energia 
de gap do GaAs, é possível determinar o valor teórico para a energia de recombinação 
elétron-buraco para o nível fundamental e para os níveis excitados do poço. A energia de 
recombinação através dos níveis fundamentais do poço é dada pela soma da energia do 
nível fundamental de elétrons (E,) e do nível fundamental de buraco pesado (E, ) com a 
energia de gap do GaAs. De modo geral podemos escrever: 


Eten) = En + En + Eg(Gaas)- (13) 


Com os resultados apresentados na tabela 2 é possível determinar o valor teórico 


para a energia de recombinação elétron-buraco para o nível fundamental e os níveis 
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excitados do poço quântico através da equação (13). 


Eceo = 0,016401 eV + 0,0021698 eV + 1,51884 eV = 1,53741 eV 
Ete1) = 0,064330 eV + 0,0021698 eV + 1,51884 eV = 1,58534 eV 
Etez) = 0,13796 eV + 0,0021698 eV + 1,51884 eV = 1,65897 eV 


Com estes resultados foi possível obter numericamente as funções de onda e a 
densidade de probabilidade nos níveis quantizados da banda de condução. 


E,=0,13796 eV 


PE )P (10º) 





“200 “00 0 100 200 





Figura 4 — (a) funções de onda (0,3 p, (x)) e (b) densidades de probabilidade calculadas para 
os 3 níveis do poço. 


Na figura 4(a) temos as funções de onda para o poço quântico, onde as mesmas 
foram multiplicadas por um fator de escala de 0,3 para uma melhor visualização destas 
individualmente. Ao analisarmos as funções de onda representadas na figura 4(a), nota-se 
que as mesmas se estendem para as regiões externas ao poço, esse é um fenômeno que 
classicamente seria impossível de ocorrer já que E < V,. Outro aspecto importante das 
regiões classicamente proibidas é que enquanto menor a energia mais rapidamente as 
funções de onda tendem a zero. 

A partir das funções de onda ao quadrado é possível obter a densidade de 
probabilidade para os níveis energéticos, mostrados na figura 4(b). Podemos notar que 
o pico que representa o nível fundamental do poço (n=0) é maior em relação aos picos 
para os níveis excitados e isso significa que a probabilidade de encontrar o elétron no seu 
estado fundamental é maior quando comparado aos estados excitados. 

Na figura 5 apresentamos o espectro de fotoluminescência (PL) da amostra 
estudada, obtido a temperatura de 8 Ke potência de 6 uW. Nesse regime de potência 
e com essa configuração experimental podemos notar a presença de um pico de maior 


intensidade denominado pico principal. No regime de potência trabalhado não foi possível 
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observar transições elétron-buraco para o primeiro e segundo nível excitado. 


T=8K Pico Principal 
Potência = 6 uiW 1,5322 eV 


Intendidade de PL (u.a.) 





1,920 | RS Pis 1,530 Rss 1,940 RE 
Energia (eV) 


Figura 5 — Espectro de fotoluminescência à temperatura de 8 K. 


A partir do espectro de fotoluminescência apresentado na figura 5 é possível obter 
a energia do pico principal que é de 1,5322 eV. Ao compararmos com o valor teórico de 
1,53741 eV temos uma boa concordância entre o valor teórico e experimental. 


5 | CONCLUSÃO 


No presente trabalho realizou-se o cálculo teórico dos níveis de energia de um 
poço de GaAs com barreiras de AlGaAs. Foram obtidos os valores teóricos das energias 
para transição elétron-buraco, possibilitando obter as funções de onda e as densidades de 
probabilidade. Comparando o valor de energia do pico principal do espectro fotoluminescência 
da amostra com as possíveis transições calculadas foi possível identificar que no regime 
de potência trabalhado não há transições de níveis excitados do poço e que a energia de 
transição teórica do nível fundamental da banda de condução para o nível fundamental 
da banda de valência está em um bom acordo com o valor experimental. Analisando as 
funções densidade de probabilidade, notamos que é mais provável encontrar o elétron no 
nível fundamental. Portanto, os resultados obtidos pelo uso de um modelo amplamente 
conhecido de mecânica quântica são muito bem exemplificados através da comparação 
com os resultados experimentais medidos por uma técnica de espectroscopia aplicada a 


um poço quântico real. 
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RESUMO: A água precipitável (AP) é uma 
importante variável atmosférica. Sua estimativa 
é de grande importância pois pode determinar 
a intensidade de um fenômeno com chuvas 
extremas ou períodos de estiagem. Em geral suas 
estimativas são feitas por várias técnicas como 
através de imagens de satélites meteorológicos, 
perfis de radiosondas e medidas por radiômetros. 
No presente estudo usamos um termômetro de 
infravermelho de baixo custo para estimar a água 
precipitável na região de Florianópolis-SC. No 
período entre Maio de 2015 e Junho de 2016 o 
termômetro de infravermelho foi usado para medir 
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INFRAVERMELHO 


a temperatura aparente da atmosfera através de 
apontamento do termômetro na direção do zênite 
celeste. As medidas foram feitas no horário de 
12:00 UTC para coincidir com as estimativas de 
Agua Precipitável (AP) feitas por radiosondagens 
do aeroporto de Florianópolis. Comparações 
entre as medidas feitas com o termômetro de 
infravermelho e com medidas por radiômetros 
mostraram forte aderência entre os dados. 
Comparações entre as medidas de temperatura 
e de estimativas de AP para dias de céu claro 
mostraram não apenas alta correlação entre 
os dados, mas permitiu também estabelecer 
uma equação empírica para estimativas de AP 
através de eventuais medidas de temperatura 
com o termômetro de infravermelho. 
PALAVRAS - CHAVE: Água precipitável, 
radiação, termômetro de infravermelho, 
Florianópolis. 


ESTIMATES OF PRECIPITABLE WATER 
FROM AN INFRARED THERMOMETER 


ABSTRACT: Precipitable water (PW) is an 
important atmospheric variable. lts estimate is of 
great importance because it can determine the 
intensity of a phenomenon such as extreme rains 
or drought periods. In general, their estimates 
are made by various techniques, such as images 
from 'meteorological satellites, radiosonde 
profiles and measurements by radiometers. In 
the present study, we used a low-cost infrared 
thermometer to estimate the precipitable water in 
the city of Florianópolis-SC region. In the period 
between May 2015 and June 2016, the infrared 
thermometer was used to measure the apparent 
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temperature of the atmosphere by pointing the thermometer towards the celestial zenith. 
The measurements were carried on at 12:00 UTC time to coincide with the estimates of PW 
made by radiosondes from the Florianópolis airport. Comparisons between measurements 
made with the infrared thermometer and measurements with radiometers showed strong 
correlation between these variables. Comparisons between temperature measuremenis and 
PW estimates for clear sky days showed not only a high correlation between the variables, but 
also allowed to establish an empirical equation for PW estimates through possible temperature 
measurements with the infrared thermometer. 

KEYWORDS: Precipitable water, radiation, infrared thermometer, Florianopolis. 


11 INTRODUÇÃO 


A água precipitável é uma variável atmosférica de grande importância. A água na 
atmosfera pode ser observada em todos os três estados termodinâmicos: vapor de água, 
àgua condensada ou partículas de gelo. As moléculas de água são um importante gás de 
efeito estufa porque podem absorver radiação térmica em várias bandas de acordo com 
seus modos de rotação e vibração (LIOU 2002). Assim, quanto maior a concentração de 
àgua mais forte é o fluxo de radiação de ondas longas da atmosfera para a superfície. 

A água precipitável (AP) pode ser estimada pela razão de mistura (gramas de 
àgua por quilograma de ar) das moléculas de água que estão presentes em uma coluna 
atmosférica. Estimativas da AP podem ser obtidas por várias maneiras, como através de 
radiosondas que são perfis medidos por balão atmosférico, radiômetros in situ e sensores 
infravermelhos. A Tabela 01 apresenta as principais técnicas usadas para medidas de AP. A 
maioria dessas medidas são caras, ou feitas esporadicamente, e não estão disponíveis na 
maioria dos locais. Portanto, encontrar maneiras de medir a AP usando sensores de baixo 


custo pode melhorar nossa capacidade de representar melhor as condições atmosféricas. 


Tipo de medida de AP Referencia 
Espectrômetro de Prisma FOWLER, 1912 


VOLZ, 1974 

(MIMS Ill 1992; BROOKS, 2001) 
(PETTIFER, 2009) 

LILJEGREN 1994 


Ocultação por GPS (GUTMAN & BENJAMIN, 2001; BEVIS et al. 
1992 


Satélite Sensor MODIS KAUFMAN & GAO, 1992 
Satélite TIROS/TOVS SODEN & LANZANTE 1995 
Radiômetro de infravermelho | (MAGHRABI & CLAY 2010) 


Termômetro de infravermelho | (MIMS et al. 2011) 





Tabela 01 — Sumario de exemplos de medidas de Agua precipitável (AP) através de vários tipos 
instrumentos. 
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Recentemente, algumas estimativas de AP foram obtidas com instrumentos de 
baixo custo, como diodos emissores de luz (LEDs) que foram usados no programa GLOBE 
(BROOKS, D. R. 2001) ou por termômetro infravermelho que foi usado para a região 
do Texas na América do Norte (MIMS et al. 2011). Usando o termômetro para estimar a 
temperatura aparente do céu no zênite, Mims et al. (2011) foram capazes de estimar a AP 
local e conseguiram produzir uma relação empírica entre a AP e a temperatura aparente 
do céu. Como moléculas atmosféricas de vapor d'água (i.e. H,O) absorvem a radiação de 
onda longa, esta radiação é re-emitida de volta para a superfície. Portanto, variações na 
concentração destas moléculas alteram o efeito estufa através de mudanças no fluxo de 
radiação que pode ser detectada por um termômetro de infravermelho. 

Neste estudo apresentamos os resultados das medições de temperatura 
atmosférica realizadas com um termômetro de baixo custo durante um ano inteiro na 
cidade de Florianópolis. Em seguida, as medidas são comparadas com dados medidos por 
radiômetros padrão e por perfis atmosféricos de radiosondagens locais. Finalmente, uma 
relação empírica é obtida para fornecer estimativas de AP de acordo com as medidas de 


temperatura. 


21 METODOLOGIA 


Um termômetro infravermelho de baixo custo (cerca de 25 dólares) foi usado para 
medir a temperatura atmosférica para o período entre maio de 2015 e junho de 2016. 
Este termômetro da marca FLUS IR-802 (Fig. 1) opera na faixa de radiação entre 8-14 
micrômetros e permite medidas de temperatura na faixa de -50 “C a 550 “CC. O termômetro 
foi usado diariamente para registrar a temperatura atmosférica aparente apontando o 
termômetro para o zênite celeste. Os registros foram realizados às 12:00 UTC (09:00 hora 
local) para coincidir com o horário de lançamento das radiosondagens feitas no aeroporto 
de Florianópolis. Além das medidas de temperatura foram anotadas também as condições 
de cobertura de nuvens no momento da medida identificando casos de céu claro ou com 


presença de nuvens. 
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Figura 01 - Ilustração de um termômetro de infravermelho (marca FLUS IR-802) usado para 
medidas da temperatura aparente da atmosfera a partir do zênite celeste. 


As medidas de temperatura atmosférica foram obtidas no campus da Universidade 
Federal de Santa Catarina (UFSC) na cidade de Florianópolis (27,5ºS; 48,5ºW) localizada 
na Ilha de Santa Catarina, no Atlântico Sul. No mesmo local do campus, foram registrados 
dados a partir de um radiômetro padrão (i.e. um pirgeômetro). Estes dados obtidos no 
LABSOLAR (Laboratório de Energia Solar) do Departamento de Engenharia Mecânica da 
UFSC foram obtidos em momento quase simultâneo para comparar e avaliar a variabilidade 
das medições do termômetro infravermelho. 

Dados atmosféricos registrados a partir das radiosondagens lançados diariamente 
por volta das 12:00 UTC (09:00 hora local) do aeroporto de Florianópolis foram usados 
para estimar a AP e comparar com a variabilidade das temperaturas registradas. Neste 
caso, a AP pode ser estimada a partir de uma integração vertical dos dados medidos em 
vários níveis atmosféricos da razão de mistura (i.e. gramas de vapor d'água por unidade 
de volume de ar seco). 

Foram obtidas curvas de tendência estatística entre as temperaturas medidas, 
os dados do radiômetro e também da AP estimada a partir dos perfis atmosféricos das 
radiosondagens. A curva de tendência obtida proporcionou uma relação empírica entre 
a temperatura medida e a AP para dias de céu claro. Essa relação estatística permite 
estimativas de AP em caso de eventual medição obtida com o termômetro de infravermelho. 


31 RESULTADOS E DISCUSSÃO 


A Figura 02 mostra uma série temporal das temperaturas registradas com o 


termômetro infravermelho. Os resultados mostram que, em geral, as temperaturas 


Física: O Universo e os Fenômenos Naturais Capítulo 2 ES 


atmosféricas aparentes são mais altas durante os meses de verão e mais frias durante 
os meses de inverno. Juntamente com todos os dados registrados, nota-se que durante 
os dias de céu claro as temperaturas são mais frias, atingindo cerca de -29 “€ em julho de 
2015. Os eventos de temperaturas mais altas que chegam a 27 “C em janeiro representam 
a presença de nuvens que aumentam a radiação descendente emitida de sua base que 
está localizada nas camadas mais baixas e mais quentes da atmosfera. 
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Figura 02 — Série temporal da temperatura atmosférica aparente (C) do céu no zênite medida 
com o termômetro infravermelho no município de Florianópolis (27,5 S; 48,5 W). 


AFigura 03 apresenta uma série temporal datemperatura registrada e sua comparação 
com a radiação quase simultânea medida com o radiômetro padrão (i.e. pirgeômetro). Os 
resultados mostram uma evolução muito próxima entre ambas as medidas. Em geral, as 
medidas de temperaturas mais altas concordam com o forte fluxo de radiação de ondas 
longas descendentes. Esses resultados mostram que o termômetro infravermelho pode 
fornecer uma boa estimativa do fluxo de radiação de onda longa descendente medida na 
superfície. 
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Figura 03 — Intercomparação entre a temperatura atmosférica (C) do termômetro infravermelho 
e a radiação (Wm?) registrada a partir de um radiômetro padrão (i.e. pirgeômetro) feita no 
LABSOLAR da UFSC. 


A Figura 04 apresenta a temperatura registrada junto com as medições de radiação e 
sua curva de tendência. Os resultados mostram uma correlação muito alta (Rº = 0,87) entre 
o termômetro infravermelho e as medições padrão do radiômetro mostrando novamente 
que o termômetro de infravermelho pode fornecer uma estimativa muito boa da condição 


atmosférica não apenas para condições de céu claro, mas também para a presença local 
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Figura 04 — Correlação entre a temperatura atmosférica (C) do termômetro infravermelho e a 
radiação (Wm*?) registrada a partir de um radiômetro padrão. 


A Figura 05 apresenta a série temporal da AP estimada a partir dos perfis das 
radiosondas e da temperatura medida com o termômetro de infravermelho. Embora ambos 
os dados apresentem a variabilidade incluindo o máximo de AP durante os meses de 


verão e o mínimo durante os meses de inverno, há vários registros discrepantes entre as 


Física: O Universo e os Fenômenos Naturais Capítulo 2 


medidas. As diferenças ocasionais ocorrem porque os lançamentos de radiosonda são 
feitos no aeroporto localizado cerca de 7 km ao sul distantes das medições de temperatura 
feitas in situ no campus da universidade. Além disso, os balões de radiosonda se movem 
devido ao vento e, portanto, registram dados longe das medições in situ e eventualmente 
poderão atingir nuvens não observadas no local dos registros de temperatura. 
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Figura 05 — Série temporal da Água Precipitada (AP) estimada a partir da radiosonda do 
aeroporto de Florianópolis e medições de temperatura infravermelha para todos os dias às 
12:00 UTC (09:00 hora local). 


A Figura 06 mostra a relação entre a AP estimada a partir das radiosondas do 
aeroporto com as temperaturas medidas com o termômetro de infravermelho para todos 
os dias. Os resultados mostram que, de fato, há uma maior correlação para temperaturas 
mais baixas que ocorrem durante os dias de céu claro e uma relação mais difusa durante 
eventos quentes nublados. 
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Figura 06 — Correlação entre a Água Precipitada (AP) estimada a partir do rádio do aeroporto 


de Florianópolis e medições de temperatura infravermelha para todos os dias. 


A Figura 07 apresenta as mesmas medidas de AP e temperatura, mas apenas para 


dias de céu claro. Os resultados mostram uma correlação muito maior (R* = 0,85) entre as 


variáveis apresentadas na Figura 06 calculada para todos os dias. 
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Figura 07 — Série temporal da Água Precipitada (AP) estimada a partir de rádiosondas do 
aeroporto de Florianópolis e medições de temperatura infravermelha apenas para dias de céu 
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A curva de tendência obtida (Fig. 08) entre a AP e a temperatura (T) apresenta a 
seguinte relação empírica: 
AP = 28.1 + 0.68T 
Esta relação empírica permite uma estimativa da AP a partir de medidas de 
temperatura com o termômetro. Portanto, este termômetro de baixo custo permite 
estimativas da AP para eventuais casos de céu claro. 
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Figura 08 - Correlação entre a Água Precipitada (AP) estimada a partir da radiosonda do 
aeroporto de Florianópolis e medições de temperatura com termômetro de infravermelho. 
Comparação com dados obtidos de MIMS et al., 2011 para região do Texas (EUA). 


4 | CONCLUSÕES E DISCUSSÃO 


Um termômetro infravermelho de baixo custo foi usado para registrar as temperaturas 
atmosféricas aparentes para um local na cidade de Florianópolis. Os resultados mostram 
que a série temporal de medições de temperatura se compara muito bem com os registros 
de um radiômetro padrão. Além disso, hã uma correlação muito alta entre as medidas de 
temperatura e a AP para os dias de céu claro. Uma comparação dessa curva de tendência 
com a obtida na região do Texas (MIMS et al., 2011) mostra valores ligeiramente mais 
baixos de AP para aquela região para condições de temperaturas mais frias. Uma vez 
que nas regiões de latitudes mais baixas (e. g. Florianópolis) a altura da tropopausa é 
mais alta, espera-se uma maior concentração de AP nessas regiões em comparação com 
as regiões de latitude mais altas (e.g. Texas). Florianópolis está localizada em uma ilha 
tropical localizada na latitude de 27,5 ºS, enquanto as medidas do Texas foram obtidas 
na latitude de 29,6 ºN. Assim, os resultados mostram que podemos esperar uma relação 


ligeiramente diferente em função da localização. Esta técnica com uso de termômetro de 
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infravermelho pode fornecer portanto, uma alternativa barata para estimativas de água 


precipitável para qualquer localidade e a qualquer instante. 
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RESUMO: Em nosso trabalho estudamos 
via Dinâmica Molecular Clássica (MD) as 
propriedades termodinâmicas em sistemas 
binários e ternários compostos por óleo leve, 
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CO, e salmoura para aplicações envolvendo 
recuperação melhorada de petróleo (RMP). 
Assumimos dois tipos de interfaces: binárias, 
onde estamos interessados na difusão das 
moléculas de CO, por uma camada variável de 
H,O ou salmoura, e atingir a fase contendo um 
modelo de óleo; e as ternárias, tal que tanto as 
moléculas de CO,, quanto a HO ou salmoura 
fazem interface com o óleo. A partir dos perfis de 
densidade para sistemas binários, foi observado 
um acúmulo de moléculas aromáticas próximo à 
interface de acordo com a espessura da camada 
de salmoura, mas nos sistemas ternários 
esse acúmulo não é observado devido à alta 
difusividade de CO, na interface direta com o óleo. 
Quando a espessura da camada de salmoura é 
maior que 4,0 nm, nossos cálculos indicam um 
aumento na tensão interfacial salmoura/óleo 
tanto com o aumento da pressão como pela 
concentração de sal. Para valores menores na 
espessura da camada de salmoura, a tensão. 
Estimamos também a diminuição na ordem de 
30% para a densidade da mistura óleo + CO, 
para sistemas ternários, quando comparada à 
densidade do óleo leve. Além disso, a salinidade 
é responsável pela estabilização da interface 
(óleo+CO,)/salmoura na formação das emulsões 
o/a. 

PALAVRAS - CHAVE: molhabilidade, interfaces, 
dinâmica molecular, recuperação melhorada de 
petróleo, emulsões. 


ABSTRACT: In our work we studied by Classical 
Molecular Dynamics (MD) the thermodynamic 
properties in binary and ternary systems 
composed of crude oil, CO2 and brine fluids for 
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the applications involving enhanced oil recovery (EOR). We assume two types of interfaces: 
binary, where we are concerned with the diffusion of the CO2 molecules by a variable layer 
of H20 or brine when achieve the oil phase; and ternary, such that both the CO2 molecules, 
as to the H20 or brine interface with the oil. From the density profiles for binary systems, an 
accumulation of aromatic molecules close to the interface was observed according to the 
thickness of the brine layer, but in the ternary systems such accumulation is not observed 
due to the high COZ diffusivity at the oil interface. When brine layer thickness is greater than 
4.0 nm, our calculations indicate an increase in the interfacial tension brine/oil with increasing 
pressure and salt concentration, but close to the interfaces containing CO2 and brine/oil 
the interfacial tension decreases with increasing pressure and salt concentration. We also 
estimate the decrease in order of 30% for the density of the mixture oil+CO2Z for ternary 
systems, when compared to the light oil density modeled in this work. In addition, salinity is 
responsible for the stabilization of the (oil+CO2)/brine interface on the formation of the w/o 
emulsions. 

KEYWORDS: Wettability, ternary interface, molecular dynamic, light oil, enhanced oil 
recovery, emulsion. 


11 INTRODUÇÃO 


A complexidade exigida nas técnicas de extração e exploração de petróleo de um 
reservatório é muito elevada e requer altos conhecimentos e estratégias em técnicas nas 
diferentes áreas da engenharia aplicada. Por considerar diversas variáveis termodinâmicas 
a serem incorporadas aos métodos atuais, as investigações em escala molecular acabam 
merecendo grande atenção [1-4]. Essas variáveis são classificadas de acordo com meios 
geológicos, pressão do reservatório de óleo, gás e ou composição da água, viscosidade 
do petróleo bruto (densidade API), migração fina, conteúdo e composição de água, 
potencial de corrosão, concentração de gás, concentração de asfalteno e parafina, entre 
outros. percebe-se que dentro dessas considerações, a presença de água é determinante 
(por exemplo, corrosão e emulsão) e geralmente ligada aos processos envolvidos na 
extração via interação de fluidos miscíveis/imiscíveis [5-9]. Sabe-se que tais interfaces 
são extremamente pequenas, da ordem do comprimento molecular, de modo que muitas 
medidas experimentais acabam sendo limitadas principalmente quando aplicadas a 
fluidos. [10-15]. Em processos que envolvem misturas de fluidos, as identificações das 
contribuições de cada componente acabam sendo gerais e não esclarecem completamente 
os mecanismos moleculares envolvidos na dinâmica de interação entre interfaces de tais 
sistemas [16]. 

O grande interesse no estudo de interfaces se reflete em técnicas de recuperação 
melhorada de petróleo (RMP) uma vez que o contato entre os fluidos é observado, onde 
grande parte do óleo contido nos reservatórios requer métodos alternativos de injeção 
de fluido [17] como uma primeira extração direta, tomando vantagem das propriedades 
umectantes entre os fluidos [18]. Sabe-se que a redução da tensão interfacial (IFT) ou 
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viscosidade do petróleo bruto por diferentes processos moleculares, envolvendo, por 
exemplo, misturas de fluidos, é o principal mecanismo dos processos RMP, uma vez que 
tais propriedades têm sua origem na interface óleo/fluido [19-21]. 

Os fluidos podem ser injetados em reservatórios de óleo por diferentes razões, 
como para manter a pressão do reservatório para melhor recuperação de óleo ou mesmo 
para tornar o óleo menos viscoso, promovendo um ganho na extração de óleo devido à 
dissolução do óleo em tais fluidos. Dentre os fluidos utilizados, a salmoura é o método 
mais comum de aplicação em reservatório devido à sua disponibilidade e manuseio. Outro 
fluido amplamente utilizado na recuperação de óleo [22,23] é o CO,, principalmente por 
seus interesses ambientais, proporcionando uma alternativa ao lançamento na atmosfera, 
por meio de sua movimentação para o reservatório de óleo [24]. Além disso, os aquíferos 
salinos podem se beneficiar do CO, armazenado no ciclo do carbono. 

Durante o processo de injeção imiscível de fluidos, o CO, é injetado na formação 
rochosa por meio da mobilização do óleo e sua movimentação para o poço de produção, 
uma vez diluído em óleo e posteriormente dilatado. [25]. Isso representa uma mobilidade 
melhorada do óleo devido à solubilidade de CO, no óleo, o que causa redução na viscosidade 
e diminuição da tensão interfacial [26]. Os principais processos imiscíveis em CO, incluem 
principalmente a injeção alternada de H,O-CO, e a injeção contínua de fluidos de CO,. A 
injeção alternada de H,O-CO, parece ser mais eficiente em aumentar a recuperação de 
óleo pesado, embora o fator de recuperação ainda seja baixo [27, 28]. Uma alternativa é o 
uso de salmoura alternada com CO, em reservatórios de petróleo e gás esgotados, onde 
normalmente há uma grande quantidade de salmoura após a inundação de H,O [29]. Isso 
leva à formação de salmoura em um sistema ternário, 

Portanto, o conhecimento do equilíbrio de fases nas misturas aquosas com 
hidrocarbonetos acaba sendo importante nos processos de recuperação de óleo. 
Conhecer a solubilidade e a volatilidade dos hidrocarbonetos é necessário para descrever 
sua distribuição de fases ao longo do processo de separação ou para auxiliar no projeto 
do equipamento de separação. Esta solubilidade mútua é amplamente dependente 
do comprimento da cadeia de hidrocarboneto e da presença de grupos polares, como 
funcionalidades aromáticas. 

Do ponto de vista da modelagem molecular em simulações computacionais, 
as misturas de H,O, óleo leve e CO, são desafiadoras devido às complexidades das 
interações entre os componentes que resultam em um comportamento muitas vezes fora 
do equilíbrio termodinâmico. Enquanto os modelos de óleo leve têm propriedades apolares, 
o CO, tem um momento de quadrupolo alto e H,O tem um grande momento de dipolo, ao 
mesmo tempo que possui ligações de hidrogênio. As misturas são, portanto, altamente 
assimétricas, tanto em termos de tamanho molecular quanto de interações de energia 
[32]. Na verdade, estudos muito limitados são observados para o sistema ternário H,O- 
CO,-óleo, especialmente o comportamento microscópico da distribuição das moléculas de 
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CO, na interface, as propriedades estruturais e dinâmicas da interface H,O/óleo são pouco 
relatadas para seus mecanismos de ação [33]. 

Diante da complexidade em se efetivar simulações computacionais em sistemas 
ternários H,O-CO,-Óleo, neste trabalho, usando simulações de dinâmica molecular 
clássica (MD), determinamos propriedades termodinâmicas como perfil de densidade, 
tensões interfaciais, carga cumulativa e coeficientes de difusão em condições típicas 
em reservatórios com o objetivo de investigar os efeitos da injeção alternada de CO, e 
salmoura em relação a sua interface com o óleo. Portanto, efeitos como a espessura da 
camada de salmoura presente entre o CO, e o óleo foram investigados pelos efeitos do 
aumento da concentração e da pressão do sal. Além disso, consideramos interfaces duplas 
(sistemas ternários) com o óleo contendo CO, e salmoura ou HO, cujo objetivo é analisar 
o comportamento da alta miscibilidade de CO, quando em conjunto com as interações de 
H,O e salmoura. 


21 METODOLOGIA 


Simulações de dinâmica molecular (MD) foram utilizadas para investigar o 
comportamento interfacial do CO,-salmoura-óleo (interfaces binárias e ternárias). 
Potenciais interatômicos precisos foram usados: moléculas de H,O foram descritas usando 
potencial SPCE / FH [34]. Os íons Na, Ca e Cl são descritas apenas por termos de van 
der Waals e Coulomb, onde esses parâmetros potenciais também são encontrados na Ref. 
[34]. Diferentes concentrações de sal foram consideradas de acordo com (NaCl/CaCl,)% 
de concentrações: (8,0/2,0), (6,4/1,6) e (4,0/1,0)%. Para moléculas de CO, foi modelado 
usando o método de ligações flexíveis [35], que é baseado no modelo EPM2 originalmente 
rígido. Para o modelo de petróleo bruto, consideramos oito tipos de hidrocarbonetos: 
alcanos [144 Hexano (HEX), 132 Heptano (HEP), 156 Octano (OCT), 180 Nonano (NON), 
96 Cicloexano (CHEX), 156 Cicloheptano (CHEP)] , moléculas aromáticas [156 tolueno 
(TOL) e 60 benzeno (BEN)]. Hidrocarbonetos foram descritos usando campos de força 
baseados em CHARMM [36]. 

As caixas de simulações para os sistemas isolados (H,O ou salmoura, CO, e modelo 
de óleo) foram configuradas como uma caixa cúbica de dimensões L x L xL eo total de 
moléculas pode ser visto na tabela 1 de acordo com as diferentes espessuras (L.) para a 
salmoura e o CO.. 
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salmoura) (nm ara CO2 (nm 
460 


4640 


Doo 
FR + 5808 
O 


O 9284 
a O I609 





2.0 
2,5 
3,0 6970 
5,0 


Tabela 1: Moléculas totais na interface binária de fluido/óleo de acordo com a espessura L do 
fluido. 


Aqui (L,xL) 


y po 
para L no modelo de célula de óleo é L, = 10,0 nm. Uma sequência de cálculos NVE 


(10,0 nm x 10,0 nm) para todos os sistemas isolados e o valor inicial 


(conjunto microcanônico), NVT (conjunto canônico) e NPT (conjunto isotérmico-isobárico) 
foram realizados para obter as densidades de equilíbrio para cada um desses fluidos em 
um conjunto de pressões e temperaturas de 150 e 300 atm, e, 300 K respectivamente. 

Neste protocolo de cálculo, para cada sistema, executamos cálculos na duração de 
1,0 ps no conjunto NVE, 10,0 ps no conjunto NVT e 20,0 ns no conjunto NPT. As condições 
de contorno periódicas foram aplicadas para todos os sistemas e interações eletrostáticas 
de longo alcance, assim como o método do espaço recíproco Particle Particle Particle- 
Mesh (PPPM) [37-39]. Em todos os cálculos, um intervalo de tempo de 0,5 fs e raio de corte 
de 12,0 À foi usado para as interações de van der Waals. A temperatura e a pressão foram 
controladas pelo termostato Nose-Hoover e pelo barostato Andersen, respectivamente. 

As interfaces binárias fluido/óleo foram construídas considerando os sistemas 
isolados equilibrados e assumindo a seguinte nomenclatura: (H,O ou Salmoura)/CO.,/óleo. 
O mesmo protocolo usado para os fluidos puros foi utilizado aqui. Após o equilíbrio, foram 
realizadas por 80,0 ns em NVT. Nesta etapa, foram calculados os perfis de densidade, 
tensão interfacial entre esses sistemas, propriedades de difusão (perfil de autodifusão e 
perfil de difusão). 

Na Figura 1 foi mostrada a configuração inicial de algumas interfaces binárias 
fluido/ óleo. Além da interface binária apresentada pelas Figuras 1(a) e 1(b), consideramos 
uma interface ternária composta por CO2 e salmoura com interface com o óleo, com a 
nomenclatura: (H,O ou Salmoura)+CO +óleo, conforme mostrado na Figura 1(c). Aqui foi 
definido como uma caixa cúbica com dimensões L x L, =5,0nmx5,0nme L =10,0 nm 
para salmoura e moléculas de COZ. Para o modelo de óleo, as dimensões definidas como 
10,0 nm x 10,0 nm x 10,0 nm. Todos os cálculos de MD foram realizados usando o pacote 
Large Atomic/Molecular Massively Parallel Simulator (LAMMPS) [40]. 
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Figura 1: Configuração inicial das interfaces binárias (H,O ou Salmoura)/CO,/óleo (após a 
fase de equilíbrio). Espessura (L ) para a salmoura e fluidos de CO,: (a) 5,0 e 10,0 nm; (b) 2,5 
e 5,0 nm, respectivamente. (c) Interfaces ternárias (H,O ou Salmoura)+CO +óleo. Todas as 
configurações em 150 atm e 300K. 


Usamos a descrição em termos de tensores de pressão para os cálculos de tensão 
interfacial usando os mesmos procedimentos vistos em [35,38]. 


1 
r=5 | de(pa(z)- pre), (1) 


—L, 


onde, as distâncias L e -L, determinam os limites da região interfacial, e p (Z) e 
P.»(Z) são OS componentes tangencial e normal do tensor de pressão, respectivamente. 

Todos os valores do coeficiente de difusão (D) apresentados aqui foram calculados 
usando a relação de Einstein, de acordo com a qual o coeficiente de difusão é calculado a 
partir do deslocamento quadrático médio: 


| N 
D= lim À 
61 N 4 


(2) 


onde r(t) é a posição desdobrada dos átomos i no tempo t, e os colchetes indicam 
uma média do conjunto, sobre todas as moléculas de soluto e todas as origens temporais. 
Também realizamos uma análise eletrostática dos sistemas, com o objetivo de 
calcular a carga iônica líquida cumulativa Q (r) nas interfaces de salmoura/((CO +óleo) de 
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acordo com, 


O(r) = [ [e.(1)-e (r)ldr, (3) 


onde c (r)ou c (r) são a distribuição de carga positiva ou negativa, respectivamente, 
definida por c (1n)=q, (1) PINRº. R é o raio máximo para o sistema, definido com a origem 
localizada no centro de massa da gota cilíndrica de salmoura, r=(Zº+y*)'* com (Zz, y) 
coordenadas cartesianas para átomos, N é o número de configurações nas simulações e | 
é o número de camadas [41]. 


31 RESULTADO E DISCUSSÕES 


3.1 Fluidos binários 


Inicialmente apresentamos resultados para as interfaces binárias representadas 
nas Figuras 1(a) e 1(b), analisando o comportamento da água ou salmoura quando se 
forma uma película ao redor do óleo evitando o contato direto CO, -Óleo. Para esses filmes, 
consideramos um estudo sistemático sobre sua variação de comprimento descrito por 
nossos modelos conforme a Tabela 1. O objetivo é verificar a permeabilidade do CO, pelo 
filme de água ou salmoura até atingir a fase oleosa e como o CO, altera as propriedades 
de molhabilidade da H,O-óleo nestes sistemas. 

Nas Figuras 2((a) a (d)), mostramos o perfil de densidade 2d ao longo da direção 
perpendicular à interface (eixo z) para interface binária de fluido/óleo a 300 K e 150 atm com 
5,0 nm de espessuras para os sistemas de salmoura. De acordo com [26], foi demonstrado 
que em uma interface imiscível fluido/óleo ocorre acúmulo de moléculas aromáticas nesta 
interface. Porém, percebe-se que no sistema que contém mais o componente CO,, ocorre 
uma dispersão de moléculas aromáticas na fase aquosa, enquanto as demais moléculas, 
predominantemente de moléculas de alcano, permanecem na fase oleosa, ver Figura 2(a) 
e 2( b). 

Além da permeabilidade das moléculas aromáticas, a permeabilidade das moléculas 
de COZ (Figura 2 (c)), não se mostra através da fase aquosa, não atingindo também a 
região oleosa. Para a densidade total de todos os fluidos, de forma a considerar também 
os efeitos da salinidade, apresentamos na Tabela 2 os valores correspondentes de 
acordo com a espessura da camada de salmoura. Em termos das concentrações máxima 
(8,0/2,0%) e mínima (4,0/1,0%) dos sais NaCl / CaCl,, observamos uma variação da 
densidade da salmoura a 150 e 300 atm de (0,99 + 0,02) g/ecmº e (1,01 + 0,02) g/em* 
respectivamente. Esses valores são valores típicos para água do mar. Para a densidade de 
CO,, considerando as variações de espessura da camada de H,O pura a 300K sob pressão 
de 150 e 300 atm o valor foi estimado em (0,85 + 0,02) g/cmº (Figura 2 (c)) e (0,92 + 0,02) 
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Figura 2: Perfil de densidade de acordo com o espessura do fluido: 5,0 nm de soluções de 
salmoura (8,0/2,0%) NaCl/CaCl,; (a) Apenas moléculas de alcano, (b) Moléculas aromáticas, 
(c) Salmoura+CO,, sem difusão de CO, na fase oleosa, e, (d) Sistema Total. Interfaces a 300K 
e 150 atm. 
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Para o nosso modelo de óleo, a densidade foi estimada em (0,76 + 0,02) g/emº e 
(0,72 + 0,02) g/cm* assumindo os valores máximo e mínimo respectivamente, da salinidade 
da salmoura, todos a 150 atm. À 300 atm, uma mudança significativa nos valores da 
densidade do óleo é observada devido à permeabilidade das moléculas aromáticas na H,O 
pura ou na fase de salmoura, (0,73 + 0,02) g/cmº e (0,70 + 0,02) g/cm*, respectivamente. 

Os resultados obtidos para os perfis de densidade de todos os sistemas binários 
foram compatíveis com os experimentais [42-45]. É possível verificar que os potenciais 
utilizados para a modelagem desses sistemas nas simulações de dinâmica molecular 
foram eficientes na captura da informação termodinâmica devido às interações entre as 


três fases (Óleo-(salmoura ou H,0)-CO,). 


Espessura ou L; do H2O Densidade (g/cm*) de acordo com a salinidade (NaCI/ 
(salmoura) CaCb%) a 300K e 150atm 
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Tabela 2: Densidades de salmoura na interface binária fluido/óleo de acordo com a espessura. 


Em relação à estabilidade da fase, ou seja, a imiscibilidade dos fluidos envolvidos, 
o aumento da concentração de sal promove uma separação de fases mais estável quando 
analisadas espessuras da camada de salmoura maiores que 3,0 nm. Considerando as 
diferentes espessuras, existe uma diferença na penetração das moléculas de CO, na 
salmoura, mas a fração no número total de moléculas de CO, é baixa, variando de 0,1 a 
1,0%. Também é observado um pequeno comportamento de acúmulo de moléculas de CO, 
na região de interface com H,O e salmoura (Figura 3), sendo mais intenso em soluções de 
H,O com espessura de 2,0 nm. 

No caso da água pura (H,0), também temos um leve acúmulo de CO, na interface. 
Para as concentrações mais altas de NaCl e CaCl,, o acúmulo de COZ não é significativo. A 
fim de investigar o acúmulo e a permeabilidade de CO, na interface com H,O ou salmoura, 
consideramos um zoom na região da interface (H,O ou salmoura)/CO, e mostrado na 
Figura 8. 
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O OS LO LS 2.0 2.5 3.0 
L (nm) 


Figura 3: Zoom no perfil de densidade de fluidos em um sistema para diferentes espessuras 
e concentrações de sal a 300K e 150 atm. Linha completa: Salmoura [H,0+(8,0/2,0%) NaCl/ 
CaCl,]; Linha pontilhada: H,O puro. Com 5,0 nm de espessuras para os sistemas de salmoura 
e H,O puro. L é a espessura da região centrada na interface. 


Para verificar os efeitos do acúmulo de moléculas nas interfaces CO(H,O ou 
salmoura)/óleo, estimamos valores da tensão interfacial salmoura/óleo de acordo com a 
espessura da camada de H,O (ou salmoura), nos sistemas de interface binária, pressão e 
concentração de sal, mostrada na Figura 4. 

Observa-se que a espessura da camada de H,O ou salmoura influencia os valores 
de tensão interfacial (H,O ou salmouray)/interface óleo. Esses resultados corroboram os 
obtidos para os perfis de densidade, pois quanto menor a espessura da camada, maior é a 
influência das interações entre as moléculas de CO, e o óleo. A explicação está relacionada 
a redução da tensão interfacial para espessura inferior a 3,0 nm, pelo fato de haver uma 
certa difusão de moléculas de CO, na fase H,O, reduzindo a tensão interfacial do (H,O ou 
salmoura)/Óleo. Supõe-se então que com o acúmulo de moléculas de CO, na interface com 
H,O ou salmoura sendo apenas um resultado das diferenças entre a tensões interfaciais, 
CO,/H,0 e óleo/H,O como visto em [46], embora os efeitos do H,O ou espessura da camada 
de salmoura, acabam eliminando essas propriedades interfaciais. 
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Figura 4: Tensão interfacial para as interfaces binárias fluido/óleo. A tensão interfacial depende 
da salinidade e das diferentes pressões a 300K. 


Quando as interfaces CO,/H,O e óÓleo/H,O estão próximas o suficiente, temos 
o acúmulo e assim resultando na diminuição da tensão interfacial, mas para interfaces 
distantes a regra é o contrário. Aqui, a tensão interfacial aumenta com o aumento da 
concentração de sal e pressão, quando a camada de salmoura possui espessura superior 
a 4,0 nm, enquanto para espessura inferior a 3,0 nm, a tensão interfacial diminui com 
aumento da concentração de sal e pressão [47-50]. Nossos resultados concordam com a 
teoria e experimental [13, 50]. 

Além de obter os valores de tensão interfacial para as interfaces bifásicas, analisamos 
a difusividade de H,O (ou salmoura) e CO, em suas fases bulk e em direção à fase oleosa 
em condições de pressão e salinidade. Em nossas simulações, foi observada a difusividade 
das moléculas de CO, em todos os sistemas estudados, a partir da suposição de diferentes 
espessuras de H,O ou salmoura entre as fases de CO, e oleosa, bem como nas interfaces 
diretas de CO, e H,O (ou salmoura) com as moléculas de óleo. 

Para quantificar a difusão das moléculas envolvidas nos sistemas apresentamos 
os resultados para os fluidos de injeção onde o coeficiente de difusão para cada ponto de 
estado foi calculado a partir de 20 simulações diferentes de 1,0 ns, ao longo da simulação 
realizada. As propriedades de transporte e os perfis de difusão ao longo da interface 
entre fluidos binários e, finalmente, comparando o comportamento desses sistemas na 
interface fluido-Óóleo de acordo com a espessura das camadas de H,O (ou salmoura). 
Assim, para compreender o efeito da interface fluido-óleo nas propriedades de transporte 
dos componentes, a variação do coeficiente de difusão (ao longo da direção perpendicular 
as interfaces fluido/óleo) foi monitorada em função da pressão e salinidade a 300K ( Figura 
5). Para todos os fluidos considerados, 

Para o sistema de interface (H,O ou salmoura)/óleo, espessuras entre 2,0 e 5,0 nm, 
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os efeitos da salinidade no perfil de difusão não foram significativos quando comparados 
ao caso H,O/óleo puro, mesmo em um sistema com pressão crescente. Com o aumento 
da concentração de sal, a solução é menos miscível com óleo, mas as propriedades de 
transporte não são afetadas em comparação com a interface H,O/óleo puro. No entanto, o 
efeito da espessura da camada de H,O que separa o CO, do óleo é relevante. Comparando 
as Figuras 5(a) e 5(c), a difusividade, pela menor espessura de H,O, é até 5 vezes maior 
do que no caso de espessura maior. 
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Figura 5: Perfil de difusão de acordo com espessura do fluido: 2,0 nm (a) H,O (linhas 
completas) e Salmoura H,0+(8,0/2,0%) NaCl/CaCl, (linhas pontilhadas); 5,0 nm (c) H,O e 
Salmoura H,0+(8,0/2,0%) NaCl/CaCl,. No sistema foram consideradas as moléculas de CO, à 
300K e 150 atm. As figuras (b) e (d) são zooms de (a) e (c), respectivamente. 


Colocando estes resultados em perspectiva, para os fluidos binários, as simulações 
sugerem uma diminuição da tensão interfacial e aumento da difusão do fluido na fase 
oleosa com o aumento da pressão quando a espessura da camada de H,O ou salmoura for 
inferior a 3,0 nm. Em suma, para esses sistemas, quanto menor a tensão interfacial, maior 
é a penetração do fluido na fase oleosa. Por outro lado, para salmoura (espessura maior 
que 4,0 nmy)/interface óleo, nossos cálculos indicam um aumento na tensão interfacial com 


o aumento da pressão e concentração de sal [51-56]. 
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3.2 Fluidos ternários 


Para as interfaces ternárias ((H,O ou salmoura)+CO ,/óleo), assumindo a configuração 
inicial da Figura 1(c), estimamos os perfis de densidade para todos os sistemas. A Figura 
6 mostram os resultados para H,O pura e salmoura (8,0/2,0% ) NaCl/CaCl,). Observamos 
imediatamente uma alta concentração de CO, na região contendo óleo (o óleo é omitido 
nas Figuras 6 (c)), devido à alta miscibilidade do CO, com o óleo. Os resultados mostram 
as mudanças nos perfis de densidade de CO, quando misturados com óleo. Comparando 
os dois casos na Figura 5, o pico de densidade de H,O puro é 0,06 g/emº maior do que para 
a salmoura na posição de 16,0 nm no eixo z da caixa computacional. 

Percebe-se também que durante a miscibilidade do CO, em óleo, tanto H,O quanto 
salmoura apresentam uma simetria tipo emulsão, devido ao contato de H,O ou salmoura 
com o novo sistema de óleo+CO, e, também, um maior número de moléculas de CO, no 
sistema contendo H,O puro quando comparado aos casos com salinidade. Na verdade, o 
último argumento é consistente com os efeitos da tensão interfacial estimada em [57,58], 
onde realizaram simulações de DM para o CO, supercritico/NaCl, CO,-Sistemas de solução 
mista de CaCl, e COZ/(NaCl + CaCl,) a 343 K e 20 MPa sob diferentes salinidades (força 
iônica). Neste caso, o efeito da salinidade demonstrou a estabilidade da solução salina 
aumentando a tensão interfacial de COZ com o discreto aumento da concentração salina. 


(a) 


d(g/em?) 


d(g/em?) 





z(nm) 
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Figura 6: Perfil de densidade para interface ternária ((H,O ou salmoura)+CO, /óleo); (a) 
Somente o óleo e (b) Salmoura (8,0/2,0%)/CO,; Configuração final para (c) H,O+CO,/óleo e (d) 
Salmoura+CO,/óleo, onde as moléculas de óleo são omitidas; (e) e (f) moléculas de óleo em 
relação a (c) e (d), respectivamente. Todas as interfaces a 300K e 150 atm. É possível observar 
uma alta difusão de CO, no óleo e H,O (ou salmoura) para todas as concentrações e pressões 
estudadas. 


Em relação à tensão interfacial CO,/óleo, procuramos explorar a influência das 
soluções salinas em reduzir ou aumentar a tensão interfacial de acordo com o processo de 
difusão. Com o erro envolvido nos cálculos da tensão interfacial na ordem de 0,1 mN/m, 
apresentamos os resultados de nossas simulações na Tabela 4. Observe uma ligeira 
variação na tensão interfacial CO,/óleo na fase ternária onde, com a presença de H,O puro 
a variação foi de 0,5/0,4 mN/m quando analisado sob pressões de 150/300 atm. Com a 
presença de salinidade, nossas simulações mostraram uma variação de 1,91,7 mN/m, 
para nossa maior concentração salina estudada. 


Fluido na Y (Este trabalho) Y (Experimental) 
interface (mN/m) (mN/m) 


de óleo 


6,4/5,7 (salmoura) * 23,0 - 1,0 [61] *** 
* - Salmoura: 8,0/2,0% NaCl/CaCl», ** - na interface decano; faixa de pressão: 1-200 atm a 300K. *** - na interface 
de petróleo bruto: faixa de pressão: 1-300 atm a 300K e 331K. 





Tabela 4: Tensão interfacial nas interfaces CO, /óleo de acordo com os sistemas de interface 
ternária a 300K e 150/300 atm 


No trabalho experimental desenvolvido em [62], foi demonstrado que em um sistema 


de salmoura-CO,, a molhabilidade diminui quando comparada com H,O-CO, puro em um 


sistema composto pela salmoura e a superfície de uma agulha, que fornecem gotículas de 
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CO, na mesma. Em outro trabalho [63], foram estudadas medições de ângulo de contato 
de equilíbrio de CO, -salmoura-óleo-sistemas de rocha mostrando aumento da tensão 
interfacial com o aumento da pressão, em temperatura constante. Isso é atribuído a uma 
maior solubilidade de CO, em uma pressão mais alta. Este ângulo de contato pode ser 
interpretado em termos de tensão interfacial CO,/óleo que neste caso se dá uma redução 
na tensão interfacial com o aumento da pressão corroborando com os cálculos teóricos. 
Percebe-se que a salinidade mostrou ser capaz de alterar levemente a tensão 
interfacial CO, /óleo na fase ternária. Podemos dizer que os materiais iônicos adsorvidos 
(sais) não causam grandes mudanças na tensão interfacial e podem aumentá-la em alguns 
casos dependendo do tipo de interface formada em um sistema ternário, mas isso faz pouco 
para facilitar a emulsificação. No entanto, sob condições apropriadas, tais sais podem 
ajudar a estabilizar o sistema, criando uma barreira eletrostática fina entre as gotículas de 
fase dispersa e alterando algumas propriedades físicas locais do sistema como a constante 


dielétrica, a viscosidade e a densidade [64]. 


% NaCI/CaCL, 


— 4.(0/1.0 
— 6.4/1.4 
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Figura 7: Carga líquida cumulativa, Q(r) com diferentes concentrações de sal a 300K e 150 
atm. 


Na verdade, ao estimar a carga iônica líquida cumulativa Q(r), a 300K e 150 atm, 
na interface da salmoura com a solução de óleo+CO, observada na Figura 7, tem-se a 
formação de uma dupla camada elétrica (EDL) na região entre 2,9 e 4,7 nm, com espessura 
de 1,8 nm, enquanto na superfície da salmoura, entre 3,6 e 4,7 nm, temos um excesso de 
carga positiva. Não houve mudanças significativas com o aumento da pressão. De acordo 
com a presença da EDL em sistemas óleo-salmoura, é visto na literatura que as emulsões 
o/a são estáveis, sugerindo que as fases óleo e H,O são distintas e imiscíveis com uma 


camada iônica na interface [65]. 
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Uma forma de caracterizar as emulsões é pelo tamanho das gotas. Como a 
emulsificação é um processo de agitação aleatória, a emulsão resultante é fornecida por 
um sistema polidisperso onde coexistem pequenas e grandes gotas. Ao observar nossas 
simulações, verificamos que a fase contendo HO ou salmoura adotou uma simetria 
cilíndrica, de modo que ao longo do eixo z observamos as circunferências de delimitação 
conforme Figura 6(c) e 6(d). A partir dessas circunferências, estimamos os diâmetros d 
ao longo de todo o eixo y centrado na circunferência de cada sistema contendo H20 ou 
salmoura. Consideramos 20 posições diferentes ao longo do eixo y e determinamos os 


valores de d de acordo com a Figura 8. 


(b) 





Figura 8: A determinação do diâmetro, d, para as gotas cilíndricas de (a) H,O pura e (b) 
salmoura (8.0/2.0%) NaCl/CaCl,. 


Os resultados das variações para d foram extremamente pequenos. Enquanto que 
para H,O puro, d = 8,70 + 0,05 nm, para a concentração de sal de (8,0/2,0%) NaCl/CaCl,, d 
= 8,80 + 0,02 nm. Porém, a simetria cilíndrica para as concentrações de sal é mais estável 
quando comparada ao HO puro. Outro fator a favor da estabilidade seria em relação à 
quantidade de CO, difundido na fase H,O ou salmoura. Em geral, o efeito da salinidade 
acaba protegendo a entrada de CO, na fase aquosa. 

Para os casos em que o CO, está em interfaces ternárias com H,O ou salmoura e 
óleo, os valores dos coeficientes de difusão de CO, na fase oleosa são semelhantes aos 
obtidos para sistemas sem H,O ou salmoura, na ordem de 16,0 x 10º m?/s a 150 atm e 20,0 
x 10º m?/s a 300 atm e não variando significativamente com a salinidade. Na verdade, todas 
as moléculas de CO, acabam entrando na fase de óleo formando uma mistura de CO, +óleo 
(ver Figura 6(e) e 6(f)) e apresentando uma diminuição na densidade do óleo e do CO,. Após 
a difusão do CO, no óleo, a densidade final do sistema CO, +óleo foi estimada em (0,58 + 
0,02) g/cmº e (0,67 + 0,02) g/cmº a 150 atm e 300 atm respectivamente, considerando nas 


estatísticas os sistemas contendo H,O puro e salmoura, observa-se uma redução de 30% 
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em média quando comparada à densidade do óleo puro. 


41 CONCLUSÕES 


As propriedades termodinâmicas das interfaces fluido-Óóleo de EOR foram 
investigadas a partir de simulações de Dinâmica Molecular sob diferentes condições de 
pressão e concentração de sal. Nós estudamos as acumulações de moléculas aromáticas 
em sistemas salmoura+CO (óleo ternário (consulte as Figuras 2 (a) e 2 (b)) e em interfaces 
binárias salmoura/CO,/óleo (Figura 2 (c)). Esse acumulo foi observado tanto para a salmoura 
quanto para a interface de CO, e diminui com o aumento da salinidade e da pressão (Figura 
3). A formação de uma interface em camadas rica em moléculas aromáticas pode ser um 
fator crucial na determinação do comportamento interfacial para tecnologias RMP. 

Para a interface salmoura/óleo, em todos os casos estudados, foi observado um 
aumento da tensão interfacial com o aumento da pressão e da concentração de sal. Na 
solução aquosa, o efeito da hidratação das moléculas de água sobre os íons, percebe-se 
um aumento da imiscibilidade entre óleo e salmoura na interface, dentro da concentração 
de sal estudada. 

Em relação as propriedades de transporte determinaram-se os perfis dos coeficientes 
de difusões na região de interface. Foi possível correlacionar os efeitos da tensão da interface 
e dos perfis de difusão na miscibilidade entre os componentes da interface. Para interfaces 
(H,O ou salmoura)/óleo não foi observada difusão da molécula de H,O dentro da fase 
oleosa. Porém, na região interfacial, observamos um discreto aumento da mobilidade com 
o aumento da pressão. Para o CO,, esses valores são significativos quando consideramos 
as espessuras (menos de 3,0 nm) das camadas de H,O ou salmoura separadas do óleo. Da 
mesma forma, com os menores valores de tensão interfacial, nossos resultados também 
sugerem uma completa miscibilidade entre os componentes acima de 150 atm para o CO.. 
Para a interface ternária entre óleo, HO e CO,, a difusão ocorre sem influência significativa 
por parte das moléculas de H,O, de modo que a densidade final da mistura de CO, +óleo é 
reduzida. Em termos de densidade do óleo puro, a redução foi de 30% em média. 

A salinidade alterou ligeiramente a tensão interfacial CO,/óleo na fase ternária, 
onde a salinidade não causa grandes mudanças na tensão interfacial, mas sob condições 
adequadas tais sais podem ajudar a estabilizar as emulsões, afetando algumas propriedades 
físicas locais do sistema, como a difusividade de moléculas e densidade. 

Portanto, a mudança na molhabilidade do óleo quando em contato com HO, salmoura 
e CO,, ocorre quando o CO, é injetado, modificando significativamente a recuperação do 
óleo. Como o CO, tende a se misturar com o óleo, ocorre um aumento de volume com 
consequente diminuição da densidade do óleo, alterando sua viscosidade e promovendo 


uma maior recuperação do óleo. 
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ABSTRACT: The increasing number of 
uncontrolled objects orbiting the Earth justifies 
efforis to observe and track them in order to 
avoid collisions among them and operational 
artificial satellites. These studies involve different 
disturbances and resonances in the orbital 
motion of these objects. Most of the cataloged 
space debris are found in low earth orbits (LEO). 
This work studies the orbital motion of space 
debris in the ISS's neighborhood. Real data 
from the 2-line element set provided by NORAD 
(North American Defense) are used to compare 
with the results obtained for orbital motions. The 
geopotential and atmospheric drag are used as 
perturbations in the SGP4 model. Figures show 
the time behavior of the orbital elements of the 
space debris in the process of reentry in the 
Earth and their distances to ISS showing that the 
collision risks of these objecis with operational 
spacecraft is growing. Solutions for the space 
debris mitigation are necessary. 
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IN LEO REGION 


KEYWORDS: Space Debris, Reentry, Collision 
Risk. 


RESUMO: O número crescente de objetos 
descontrolados orbitando a Terra justifica os 
esforços para observá-los e rastreá-los, a fim 
de evitar colisões entre eles e satélites artificais 
operacionais. Estes estudos envolvem diferentes 
perturbações e ressonâncias no movimento 
orbital desses objetos. A maior parte dos detritos 
espaciais catalogados é encontrada em órbitas 
baixas da Terra (LEO). Este trabalho estuda 
o movimento orbital de detritos espaciais na 
vizinhança da ISS. Dados reais do grupo de 
elementos 2-linhas fornecidos pelo NORAD 
(North American Defense) são usados para 
comparar com os resultados obtidos para os 
movimentos orbitais. O geopotencial e o arrasto 
atmosférico são usados como perturbações no 
modelo SGP4. Figuras mostram o comportamento 
no tempo dos elementos orbitais dos detritos 
espaciais no processo de reentrada na Terra e 
as distâncias destes para a ISS mostrando que o 
risco de colisão destes objetos com espaçonave 
operacional está crescendo. Soluções para a 
diminuição de detritos espaciais são necessárias. 
PALAVRAS  - CHAVE: Detritos Espaciais, 
Reentrada, Risco de Colisão. 


11 INTRODUCTION 


Since the launch of the first satellite, 
several space missions have been realized 
putting objects around the Earth. In this way, 
several studies are important to preserve the 
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operability and integrality of the operational artificial satellites, considering the increasing 
number of distinct objects in the space environment offering collision risks and, consequently, 
possible lost of mission (OSIANDER and OSTDIEK, 2009). 

The International Space Station (ISS) is a great opportunity to use a research 
platform in space. An international parinership of space agencies of the United States 
(NASA), Russia (FKA), Japan (JAXA), Europe (ESA), Canada (CSA), among other, provides 
the operation of the ISS with astronauts since 2000 (ISS, 2015). The ISS is in Low Earth 
Orbits, in the same region of most of the space debris orbiting the planet. In this way, 
several studies are important to preserve the operability and integrality of the space station 
and operational artificial satellites, considering the increasing number of distinct objects in 
the space environment offering collision risks and, consequently, possible lost of mission 
(KLINKRAD, 2006; MECHISHNEK, 1995). 

The constant monitoring of the ISS's orbit is needed, as well as the monitoring of 
the orbits of the more than three thousand Earth's satellites currently in operation is also 
required. The monitoring of temporal variation of the orbital elements of space debris in the 
LEO's region, as proposed in this paper by using of SGP4 model, can be useful of planning 
evasive maneuvers to avoid collisions and a premature reentry of spacecrafts (including the 
ISS). Another example of situation to be avoided is the collision between space debris and 
satellites that have reached the end of their useful life and are in a forced or planned reentry 
maneuvers for they do not become space debris too. In this case, collisions would generate 
more debris into LEO's region yet. Analyses of temporal variation of the orbital elements of 
space debris can also avoid these unwanted events. 

The orbital dynamics of cataloged objects can be analyzed using the 2-line element 
group, or TLE of the NORAD (North American Defense). The TLE data are composed of 
seven parameters and epoch (HOOTS and ROEHRICH, 1980; SPACE TRACK, 2019). 

Figure 1 shows that most of objects are in the region 13 <n (rev/day) < 15 by the 
histogram of the mean motion of artificial satellites and space debris in LEO region. 

Synchronous satellites in circular or elliptical orbits have been extensively studied 
in the literature, due to the study of resonant orbits characterizing the dynamics of these 
objects (SAMPAIO et al., 2014; 2016). 

In this work, a study is done about objects with orbital motions around the Earth in 
LEO region. These space debris reentry in the Planet in some moment. The TLE data of the 
cataloged objects are used and the orbit propagation is done considering the effects of the 
Geopotential and Atmospheric Drag. The orbital motions of space debris are studied in the 
neighborhood of the ISS - International Space Station. The results show that the collision 
risk of space debris with the ISS is high and solutions to avoid these events and for space 
debris mitigation are necessary. 
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Figure 1: Histogram of the mean motion of the cataloged objects in LEO region. 


In the next section, a propagation of the orbit is presented showing solutions for the 


time behavior of the orbital elements. 


21 PROPAGATION OF THE ORBIT 


In this section, a method used to propagate the orbit is described. The initial data of 
the space debris are taken from the real orbital motions shown in the Two-Line Elements 
described before (see Figure 1). First, the real data are corrected and then propagated 
in time using the SGP4 model. More details about this method is shown in (HOOTS and 
ROEHRICH, 1980). The SGP model was used by the National Space Surveillance Control 
Center and after by the United States Space Surveillance Network with the goal to track the 
objecis in orbital motions. 

A previous study, described in Sampaio et al. 2014, shows a large research in 
resonant regions with several objects, including space debris. Results show that most of 
objects are in the region 13 < n< 16, where n is the mean motion. The International Space 
Station is in the related region and studies about proximities with space debris, safety orbits 
and collision risks are natural consequences of the development of the work. 

The propagation method shown in the present work considers effects of the 
geopotential and atmospheric drag, see Eqgs. 2.7 to 2.19, while previous works, see Sampaio 
et al. 2014 and Sampaio et al. 2016, they consider only corrections in the TLE data and 
effecis of the geopotential in the orbital dynamics of space debris shown by Egs. 2.1 to 2.6. 

The corrections in TLE data shown by Egs. 2.1 to 2.6 are necessary because the 
Earth's atmosphere cause differences in the data produced by the sensors. Considering n, 


the mean motion of the 2-line, the semi-major axis a, is calculated by Eq. (2.1), 
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m=(Vu/ny)S, (2.1) 


where // is the Earth gravitational parameter, /=3.986009 x 10'* mº/s*. Using a,, the 
parameter O, is calculated by the Eq. (2.2), 


aé (3cos*(1)-1) 


3 
d1=4)2 3 “eo (2.2) 


where a, is the Earth mean equatorial radius, a =63/8.135 km, J, is the second zonal 
harmonic, J,=1,0826 x 109, e is the eccentricity and | is the inclination of the orbital plane 
with the equator. 

Now, the new semi-major axis a, is defined using 5, from the Eg. (2.2) (HOOTS; 
ROECHRICH, 1980), 


ao=01[1-56, — 8," 6,7], (23) 


81 


and a new mean motion nº and semi-major axis a” are found by the parameter 5,, 


3, aí (3cos“(1)-1) 








2 (ue? (2.4) 
Po MM Mo. 

no = E 2.5) 

aí, = — (2.6) 


The secular effects of gravitation and atmospheric drag are included in the next 
equations. 
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3Jpag(-1+302) | 3JZag(13-7802+13780?) 





M'=M+[1+ mar ri nó(t— to), (2.7) 

ro 3J,2a2(1-5602) | 3)Zaé(7-11402+3950*)  15],ag(3-3602+490*)] , 
O=wm+ [O No(t — to) (2.8) 

ro 3Jpa20 | 3J2aé(40-1903)  15],ag0(3-702)] , 

Q =0+[—-S + o — no(t — to) (2.9) 
Sw! = B'Ca(cos(w))(t— to) (2.10) 

, * e ' 3 3 
8M' = —2 (go — JB ot=l( + ncos(M))* — (1 + ncos(M)) | (2.11) 
Mº=M'+ 6w' + ô6M' (2.12) 
wo" =w — Sw — 6M' (2.13) 

PO 21 noJzaêô e 2 

O" = Q 7 apê C(t—to) (2.14) 
e” =e-B'Ci(t—to) — B'Cs(sin(M") — sin(M)) (2.15) 


a” =ao[1-Ci(t—to)-D>(t— to) — Ds(t— to) —-Da(t-to)] (2.16) 


K=M+o"+Q" +mns Ecs(t — to)? + (D, + 2C))(t— to)? + =(3D; + 12C,D, + 


10CF)(t— to)! + = (3Dy + 12C,D3 + 6D3 + 30CZDo + 15C)(t— to)*). (2.17) 


Br=V(1-e2, (2.18) 


n'= E (2.19) 


where w is the argument of pericentre, OQ is the longitude of the ascending node and 
M is the mean anomaly. (tt ) is the time since epoch, B' is the drag coefficient and J, the 
fourth gravitational zonal harmonic of the Earth. 
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Observe that Eqgs. 2.7 to 2.19 show transformations, (a, e w, O, M) -> (a”, e”, 
w”, O”, M”?, in the classical orbital elements and in other variables, as mean motion n. 
These transformations represent corrections in TLE data and effects of geopotential and 
atmospheric drag are included in orbital motions of space debris. 

Some modifications are considered in the equations and used for objecis in the 
process of reentering the Earth. For values for the perigee between 98 km and 156 km, the 
value of the constant is 1.01222928. So, itis possible to rewrite, 


e+=ao(l-e)-€E+a. (2.20) 


The value of e is changed again, when the perigee is below 98 km, 


esx=+a, (2.21) 


where X=Earth radii / kilometers = 6378.135. 


Considering the changes in the constant, the term (q -e)* is replaced by 
et — 2 AMA Len + 
(go — e+)* = [[(go — JJ” + e— ex] (2.22) 


where q, is a parameter for the SGP4 density function. 

The terms used in Egs. 2.7 to 2.17 0,0, B, mn, C2, C1, C3, Ca, Cs, D2, D3 and D, 
are described in Hoots and Roerich (1980), with the appropriate values of and (q -e)”. 

In the next section, the orbital motion of the International Space Station, ISS, is 
studied. 


31 ORBITAL MOTION OF THE ISS 


Atheory to study orbital motion of objects in LEO region can be introduced considering 
perturbations caused, by the Geopotential and atmospheric drag. The orbital dynamics of 
the International Space Station is propagated using the equations presented in the second 
section. 

Figures 2 to 5 show the time evolution of the classical keplerian elements of the ISS 
in 2018. 

The orbital dynamics of the ISS is constantly monitored and corrections can be done 
to avoid the process of reentry in the Earth or collisions with space debris. 


In Figure 2, itis possible to observe that in some time the ISS begins the process 
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of reentry needing corrections in your orbit. Several objects are in the same region of the 
ISS in Low Earth Orbits and orbital maneuvers are necessary to avoid collisions with space 
debris and the process of reentry in the Planet. 

In the next section, the orbital motions of space debris in process of reentry are 
analyzed. 
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Figure 2: Orbital motion of the ISS corresponding to January/2018 to December/2018: Time 
behavior of the semi-major axis. 
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Figure 3: Orbital motion of the ISS corresponding to January/2018 to December/2018: Time 
behavior of the eccentricity. 
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Figure 4: Orbital motion of the ISS corresponding to January/2018 to December/2018: Time 
behavior of the inclination. 
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Figure 5: Orbital motion of the ISS corresponding to January/2018 to December/2018: Time 
behavior of the argument of pericentre. 


41 REENTRY IN THE EARTH 
In this section, the orbital motions of space debris are analyzed. The objecis, in the 
process of reentry, are faster when near the Earth's surface. The space debris FENGYUN 


30158 has a prediction of reentry for May/2018. 

Figures 6 a) and 6 b) show the time behavior of semimajor axis and the time behavior 
of eccentricity, respectively, of space debris Fengyun 30158. In some months, the process 
of reentry is completed, when the object reaching the Earth's surface, considering the 
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Earth's radius. 

Figures 7 a) and 7 b) show the time behavior of semimajor axis and the time behavior 
of eccentricity, respectively, of space debris Techedsat 6 43026. In some months, the 
process of reentry is completed in May 22, 2018. 

Table 1 shows some objects in process of reentry. The space debris analyzed can 
be fragments, tools, rocket bodies from space missions ISS, IRIDIUM 33 and DELTA 1. The 
predictions are for 2018 and 2019 and they can represent a safety planning and avoid risks 
to the population that collisions with the Earth can generate. 


Cataloged number | Reentry prediction 
ISS deb 43206 November 16, 2018 


IRIDIUM 33 deb May 29, 2018 
DELTA 1 deb September 04, 2019 





Table 1: Space debris in process of reentry. 
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Figure 6: Orbital motion of space debris FENGYUN 30158 corresponding to January/01/2018 to 
May/18/2018: a) Time evolution of the semimajor axis and b) Time evolution of the eccentricity. 
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Figure 7: Orbital motion of space debris Techedsat 6 43026 corresponding to January/01/2018 
to May/22/2018: a) Time evolution of the semimajor axis and b) Time evolution of the 
eccentricity. 


In the next section, the orbital motions of some space debris are compared with the 
orbital dynamics of the ISS verifying possible collision risks. 


5 | RANGE AND COLLISION RISK 


In this section, the orbital motions of space debris and the International Space 
Station are compared, analyzing the moments of proximities between the orbits. These 
space debris are in the neighborhood of ISS and can offer collision risks. An impact in the 
ISS can damage the spaceship in several levels, depending on the size of the object that 
collided and the relative velocity between them. 

The objects, in the process of reentry, are faster when near the Earth's surface. 
The space debris Sentinel 41803 has a prediction of reentry for September/2017, showing 
more situations of proximities with the ISS; while SL-18 28961 and Pegasus 24457 had 
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predictions to reenter in January/2018 and March/2018, respectively, what was confirmed 
after by real data. Figure 8 shows the orbital motions of space debris Sentinel 41803 and 
SL-18 28961 crossing the orbital motions of the ISS, by the time behavior of the semimajor 
axis. 

Distances between the ISS and space debris Sentinel 41803 can be observed in Fig. 
9, space debris SL-18 28961 in Fig. 10 and involving the ISS and Pegasus 24547, in Fig. 11. 
Both objects starts the orbital dynamics in the neighborhoods of ISS, but, after some time, 
the process of reentry is completed for objects Sentinel 41803 and Sentinel 41803. See Fig. 
8 where can be observed, probably, the objects reaching the Earth's surface, considering 
the Earth's radius. 

Analyzing the time behavior of the distances between space debris and operational 


artificial satellites is possible to know about collision risks and if is necessary some orbital 
maneuver. 
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Figure 8: Time behavior of the semimajor axis of the ISS (ao=6782.44 km) and space debris 
SL-18 28961 (ao=6836.26 km) and Sentinel 41803 (ao=6960.46 km). 
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Figure 9: Time behavior of the distance between the orbital motions of the ISS and space debris 
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Figure 10: Time behavior of the distance between the orbital motions of the ISS and space 
debris SL-18 28961. 
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Figure 11: Time behavior of the distance between the orbital motions of the ISS and space 
debris Pegasus 24547. 


61 CONCLUSIONS 


In this work, the orbital dynamics of space debris are studied. There are several 
operational artificial satellites in Low Earth Orbits, in the same region of most of the space 
debris orbiting the planet. So, several studies are important to preserve the operability 
of the space station, considering the increasing number of distinct objects in the space 
environment offering collision risks and, consequently, possible lost of mission. 

The orbital motion of objects are propagated, considering perturbations of geopotential 
and atmospheric drag. Real data from the Two Line Elements Set of the NORAD are used 
as initial conditions to the time evolution of the orbital elements. The orbital dynamics of the 
ISS is compared with space debris and several moments of proximities are indicated. 

The process of reentry of Fengyun 30158 and Techedsat 6 43026 in the Earth are 
observed when the time behavior of the semi-major axis of the space debris are studied. 

Distances between the orbital motion of the ISS and the space debris can provide 
more security and informations about collision risks. 

One can observe that the collision risk of space debris with the International Space 
Station is growing and solutions for the space debris mitigation are necessary. 


As a future work, strategies to reduce space debris in LEO region can be studied. 
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RESUMO: No ensino básico e superior, os 
professores ao trabalhar as propriedades 
elásticas dos corpos, usam como exemplos 
o comportamento linear das molas. Esse 
comportamento é explicado através da 
denominada “Lei de Hooke”, consistindo que 
a força de deformação aplicada à mola é 
diretamente proporcional à elongação produzida 
na mola. No entanto, os modelos, apresentados 
em alguns livros didáticos, da associação em 
série e em paralelo de duas molas, predizem 
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a constante elástica equivalente como a soma 
das inversas das constantes elásticas, e a soma 
das constantes elásticas das molas individuais, 
respectivamente. Esses modelos carecem de 
adequada discussão sobre sua construção e 
não são testados empiricamente. No caso da 
associação em paralelo, a hipótese assumida no 
modelo tradicional é a igualdade das elongações 
de ambas as molas, o que não se ajusta com 
a realidade. No presente artigo, propomos 
uma expansão desse modelo tradicional para 
determinar a constante elástica equivalente de 
duas molas associadas em paralelo, na qual é 
modificada sua hipótese básica sob um enfoque 
energético. Para isto, seguiremos a abordagem 
da modelagem científica de Bunge, adotando 
um objeto-modelo e o incorporando na teoria da 
Mecânica Newtoniana. Logo, o modelo teórico 
proposto é validado empiricamente. 

PALAVRAS - CHAVE: Associação de Molas. 
Energia Potencial Elástica. Lei de Hooke. 
Modelagem Científica. 


SCIENTIFIC MODELING: 
CONSTRUCTION OF THE THEORETICAL 
MODEL OF TWO ASSOCIATED SPRINGS 

IN SERIAL AND PARALLEL UNDER THE 
ENERGY APPROACH 


ABSTRACT: In basic and higher education, 
teachers, when working on the elastic properties 
of bodies, use the linear behavior of springs as 
examples. This behavior is explained through 
the so-called “Hooke's Law”, consisting that 
the deformation force applied to the spring is 
directly proportional to the elongation produced 
in the spring. However, the models, presented in 
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some textbooks, of the association in series and in parallel of two springs lack adequate 
discussion about their construction and are not tested empirically. In this article, we propose 
an expansion of the traditional model to determine the equivalent elastic constant of two 
associated springs in parallel, in which its basic hypothesis is modified under an energetic 
approach. For this, we will follow the approach of Bunge's scientific modeling, adopting a 
model object and incorporating it into Newtonian Mechanics theory. Therefore, the proposed 
theoretical model is validated empirically. 

KEYWORDS: Springs Association. Potential Elastic Energy. Hooke's law. Scientific Modeling. 


11 INTRODUÇÃO 


No reconhecido dicionário Aurélio (2004), existem diversas definições da palavra 
modelo, dentre as quais, as mais próximas na área das ciências exatas são: “conjunto 
de hipóteses sobre a estrutura ou o comportamento de um sistema físico pelo qual se 
procuram explicar ou prever, dentro de uma teoria científica, as propriedades do sistema” 
(Ibid., p. 1345), ou, “Representação simplificada e abstrata de fenômenos ou situação 
concreta, e que serve de referência para a observação, estudo ou análise” (Ibid., p. 1345). 
Para esclarecer, detalhar e ampliar esses conceitos de modelo nas ciências, seguiremos 
a abordagem da modelagem científica segundo o físico-matemático e filósofo da ciência 
Mario Bunge (MACHADO; BRAGA, 2019). Ele define (BUNGE, 1974, p. 16): “Um modelo 
teórico é um sistema hipotético-dedutivo que concerne a um objeto-modelo, que é, por sua 
vez, uma representação conceitual esquemática de uma coisa ou de uma situação real ou 
suposta como tal”. 

Segundo Bunge (1974), existe uma variedade de escolhas de objetos-modelo m 
que resultam em modelos teóricos (também denominadas como teorias específicas Ts) 
diferentes de um mesmo sistema, sendo o modelo mais simples a caixa negra dotada 
somente de entrada e saída, que constituirá uma representação do funcionamento global 
do sistema em estudo. Ou uma caixa mais elaborada como a caixa cinzenta ou translúcida 
que além da entrada e saída é munida de mecanismo interno com estados internos. 


Segundo essa escolha, Bunge explicita (1974, p. 21): 


A escolha entre estes diversos objetos-modelo e os modelos teóricos 
correspondentes dependerá do objetivo do investigador. Se se trata apenas 
de manejar um sistema, então uma caixa negra poderá bastar; mas se se 
quer compreender o seu funcionamento, seja porque se quer dominá-lo ou 
modificá-lo então não se poderá deixar de imaginar modelos mais ou menos 
profundos, gozando de apoio de teorias gerais bem como do apoio de 
experiências novas (lbid., p. 21). 


Em relação às noções de objeto-modelo e modelo teórico na ciência fatual Bunge 
faz uns esclarecimentos (1974, p. 31): 
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Preocupar-nos-emos com objetos-modelo e modelos teóricos como esboços 
hipotéticos de coisas e fatos supostamente reais. Assim um fluido pode ser 
modelado como um contínuo dotado de certas propriedades e semelhante 
objeto-modelo pode ser enxertado em uma das várias teorias gerais, digamos 
a mecânica clássica ou mecânica relativista geral. Do mesmo modo é possível 
modelar um organismo de aprendizagem como uma caixa negra equipada 
com determinados terminais de entrada e saída e pode-se desenvolver este 
objeto modelo em um sistema dedutivo hipotético. Em qualquer dos casos 
produz-se uma teoria específica ou modelo teórico de um objeto concreto. O 
que se pode submeter a provas empíricas são tais modelos teóricos: as teorias 
gerais despreocupadas com particulares permanecem incomprováveis, 
a menos que sejam enriquecidas com modelos de seus referentes. E os 
objetos-modelo mantêm-se estéreis a não ser que sejam introduzidos ou 
desenvolvidos em alguma teoria (Ibid., p. 31). 


A seguir reproduzimos no quadro 1, dos autores Brandão et al. (2011), na qual 
ilustra com alguns exemplos da Física o processo de modelagem científica segundo Bunge 
(1974). Ressaltando que a primeira coluna denominada de situação a ser modelada é o que 


denominamos objeto - concreto a ser estudado. 


e rag Objeto-modelo Modelo teórico 
. : Modelo de fluido 
Escoamento da água no Fluido continuo sem viscosidade ideal 


interior de uma Mecânica dos Fludos 


tubulação Fludo continuo com viscosdade ata ua 
viscoso 


; É Sistema termicamente isolado de Mecânica Estatistica e | Modelo de gás 
Certa quantidade de gás ideal cléesico 
contida num recipiente 


fechado per-eitamente elásticas com as paredes | Mecânica Estatística e | Modelo de gãs 
E Mecânica Quântica | ideal quântico 


partículas que interagem via colisões Mecânica Clássica 


Comportamento da 
materia em nível 
microscópico 


Mecânica Clássica e Modelo atômico 
Eletromagnetismo de Rutherford 
Sistema de particulas sujeitas a uma 
força central. Modelo 
Mecânica Clássica gravitacional de 
Newton 


Movimento dos planetas 
do Sistema Solar 





Quadro 1: Exemplos de situações do processo da modelagem científica na Física 


Observa-se, que um mesmo objeto-concreto, escoamento da água no interior de 
uma tubulação, pode ter mais de um objeto-modelo que ao ser inserido na teoria geral da 
Mecânica dos Fluídos resulta em dois modelos teóricos. Por sua vez, um objeto-concreto, 
como o gás contido num recipiente fechado, representado por um objeto-modelo, pode ser 
enxertado em duas teorias gerais resultando em dois modelos teóricos. 

Dentro deste contexto, elasticidade é uma propriedade que ocorre em alguns 
materiais quando são deformados por uma força de compressão, de tração ou flexão, e 
voltam ao seu estado inicial quando suprimida essa força deformadora. No ensino de Física, 
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em geral, é usado molas no estudo das propriedades elásticas a partir da Lei de Hooke, que 
relaciona a força restauradora da mola com a sua deformação, dada por (NUSSENZVEIG, 
2013): 


P=-kã, (1) 


sendo, k a constante elástica da mola e X é a deformação ou a elongação da mola. 

No geral, os livros de ensino de Física ao se referirem do termo “associações”, O 
usam em situações variadas como, resistores, capacitores e molas, seja em série ou em 
paralelo. Nos quais se constroem modelos simplificados, sem especificar detalhes, com o 
auxílio de analogias entre essas associações que algumas vezes resultam inadequadas, 
com hipóteses e aproximações pouco explícitas, e sem testes empíricos. Por exemplo, 
no modelo da associação em paralelo de duas molas, alguns livros didáticos apresentam 
que a constante elástica equivalente é dada como a soma das constantes elásticas das 
molas individuais (ATAM, 1998; PARANA, 1998; YOUNG; FREEDMAN, 2008). A predição 
desse modelo se ampara na hipótese de igualdade de elongações de ambas as molas, e 
é desprovido de testes empíricos. Em recente trabalho, Campomanes, Heidemann, e Veit 
(2016) mostraram que esse modelo tradicional é válido somente quando as constantes 
elásticas de ambas as molas são iguais. No entanto, a hipótese de igualdade de elongações 
malogra quando as molas individuais apresentam constantes elásticas diferentes. Logo, é 
necessária a expansão do modelo tradicional. 

Neste trabalho, estudaremos a associação de duas molas associadas em série e em 
paralelo. Para as molas associadas em série faremos a contrastação empírica do modelo 
tradicional. Para o caso das molas associadas em paralelo, para constantes elásticas 
diferentes, é proposto um modelo teórico na qual como consequência de sua hipótese 
resulta uma constante elástica equivalente diferente do modelo tradicional. Para ambas as 
situações, seguiremos a abordagem da modelagem científica de Bunge (BUNGE, 1974; 
MACHADO; BRAGA, 2019), considerando um objeto-modelo que represente o objeto 
concreto (real), idealizando este através de simplificações e postulando uma hipótese sob 
enfoque energético e, em seguida, o objeto-modelo é inserindo à mecânica clássica como 
teoria geral. Posteriormente, este modelo teórico é testado empiricamente. 


2 | CONSTRUÇÃO DO MODELO TEÓRICO DE DUAS MOLAS ASSOCIADAS 
EM SERIE 

De acordo com Bunge (1974), o objeto-modelo é uma representação conceitual 
simplificada de seu referente, o objeto-concreto. No presente caso, o objeto-concreto 
é representado pela Figura 1, em que duas molas, de constantes elásticas k, e k, são 
colocadas em serie com uma carga pendurada no extremo livre da associação de molas. 
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Figura 1: Duas molas com constantes elásticas k, e k, associadas em serie, com uma carga 
aplicada na extremidade livre da associação. 


A Figura 2.b representa o objeto-modelo que considera algumas características do 
objeto-concreto representado na Figura 1, onde são feitas aproximações, simplificações e 
hipóteses. Considera-se que: 


) As massas das molas, da haste, e as pequenas vibrações das molas são 
desprezíveis; 


li) A associação de molas em série é representada por uma mola denominada 
equivalente de constante elástica k, (ver Figura 2.c). 


Iii) À elongação da mola equivalente “x” é dada pelo deslocamento entre a posição 
inicial do extremo livre da associação sem carga (posição x=0, linha horizontal 
tracejada da figura 2) e o ponto de suspensão da carga pendurada no extremo livre 
da associação quando as molas estão deformadas (seta central da Figura 2.b); 


iv) Os comportamentos das molas 1, 2 e equivalente obedecem a Lei de Hooke, ou 
seja, F=kx, F,=kx,e F =F.=kx,. Sendo F, a força aplicada à associação de 
duas molas, a mesma à aplicada à mola equivalente. 


v) A soma da energia potencial elástica das molas 1 e 2 é igual a energia potencial 
elástica da mola equivalente (diferentemente da soma das elongações, x, +X, = X,, 
dos livros didáticos). 


Logo; 


Física: O Universo e os Fenômenos Naturais Capítulo 5 2 o. 


= hyx7 + >kox2 =>Kkçx£ (2) 


Ki K 


K> ko 


Figura 2: (a) Representação de duas molas associadas em série sem carga, (b) representação 
do objeto-modelo proposto a partir do objeto-concreto (Figura 1) com carga aplicada F no 


extremo livre da associação de molas, (c) mola equivalente que representa as duas molas com 
a mesma carga F e elongação à representada pela parte (b) desta Figura. 


A equação (2) é a hipótese principal do presente modelo. Isolando as elongações 
da hipótese (iv), temos: 


a nt ot 2 PE 
o nm (3) 


Substituindo Eq.(3) na Eq.(2): 


que É 


k1 ks kE 


Por outro lado, a força deformadora F= F.. é a mesma força que age na mola 1 e na 
mola 2, ou seja: F = F,=F.=Fe o substituindo na Eq. (3.a), temos: 
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=+> (4) 


A Eq.(4) é a predição do valor da constante elástica equivalente do modelo teórico. 


2.1 Contrastação do Modelo Teórico 


Foram coletados os dados de Força e elongação para a carga pendurada de acordo 
a configuração experimental representada na Figura 1, como também das molas de 
maneira separada, a incerteza na medida da elongação equivalente, considerando o efeito 
de paralaxe, foi de 0,002 metros. Posteriormente os dados são plotados como mostra a 
Figura 8. 


A mola 1k, 
O mola 2k, 
* mola 1+2 ks 


k,= (8,2 + 0,4) N/m 
k= (13,4 + 0,2) Nim 
kexp= (5,3 £ 0,5) Nim 
ke= (5,1 + 0,4) N/m 





0,00 0,05 0,10 0,15 0,20 0,25 
Elongação (m) 


Figura 3: Força aplicada em função da elongação para as molas 1, 2 e equivalente. Nesta 
Última, as molas 1 e 2 são associadas em série. 


Nessa figura observa-se as curvas de ajuste linear para a determinação da constante 
elástica da mola 1, mola 2 e das molas associadas em série. Sendo esses valores de k, = 
(8,2 + 0,4) N/m, k,= (13,4 + 0,2) N/m e Koxo = (5,3 + 0,5) N/m, respectivamente. 

Enquanto que, o valor predito pela Eq. (4), deduzida dentro do modelo teórico, é k, 
= (5,1 + 0,4) N/m. Portanto, considerando a incerteza da predição, podemos inferir que o 


resultado experimental corrobora empiricamente o modelo teórico. 
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81 CONSTRUÇÃO DO MODELO TEÓRICO DE DUAS MOLAS ASSOCIADAS 
EM PARALELO 

No presente caso, o objeto-concreto é representado pela Figura 4 (Figuras 4.b e 
4.c), em que duas molas de constates elásticas diferentes (k, z K,), inicialmente de mesmo 
comprimento como mostra a Figura 4.a, são associadas em paralelo com uma haste as 
unindo nas suas extremidades inferiores com uma carga pendurado num ponto da haste. 





(a) (b) (c) 


Figura 4: (a) Duas molas com constantes elásticas diferentes (k,  K,), de igual comprimento 
associadas em paralelo com uma haste de alumínio unindo suas extremidades; (b) com uma 
carga aplicada no centro da haste; e (c) a carga aplicada próxima da mola 1. 


A Figura 5.a representa o objeto-modelo que considera algumas características do 
objeto-concreto representado na Figura 4.c, onde são feitas aproximações e hipóteses: 
i) As massas das molas, da haste, e as pequenas vibrações das molas são 
desprezíveis. 


li) As molas 1 e 2 estão na posição vertical e ambas são paralelas. 


Iii) A associação de duas molas é representada por uma mola denominada 
equivalente de constante elástica k, (ver Figura 5.b). 


Iv) a carga aplicada na haste não é necessariamente aplicada no centro deste (ver 
Figura 5.a), sendo L, a distância entre o extremo inferior da mola 1 e o ponto de 
aplicação da carga (L = L +L,, sendo Lo comprimento da haste). 
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v) À elongação da mola equivalente XE é dada pelo deslocamento entre a posição 
inicial da haste sem carga (posição x=0, linha horizontal tracejada cor vermelha) 
e o ponto de suspensão da carga pendurada na haste quando as molas estão 
deformadas (seta central da Figura 5.a). 


vi) Os comportamentos das molas 1, 2 e equivalente obedecem a Lei de Hooke, ou 
seja, F=kx, F,=kx,e F=F.=kx,. Sendo F, a força aplicada à associação de 
duas molas, considerada a mesma aplicada à mola equivalente. 

vii) Assume-se que a soma da energia potencial elástica das molas 1 e 2 é igual 


a energia potencial elástica da mola equivalente (diferentemente da hipótese de 
igualdade de elongações, dos livros didáticos). 


Logo: 


=hyxi + kox2 = —kpxf (5) 






—“”Mola 
Equivalente 


É 
(b) 


Figura 5: (a) Representação do objeto-modelo proposto a partir do objeto-concreto (Figura 4.c) 

de duas molas associadas em paralelo com elongações diferentes x,x x2, com carga aplicada 

F em qualquer ponto da haste, (b) mola equivalente que representa às duas molas da parte (a) 
com a mesma carga F e elongação x. 


Ao igual da associação de molas em série, a equação (5) é a hipótese do presente 
modelo proposto. Logo, enxertamos o objeto-modelo numa teoria geral, neste caso a 
mecânica Newtoniana. 


Fazendo o Diagrama de Corpo Livre (DCL) da haste (ver Figura 6). 
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Ly 
L> 


F 
Figura 6: DCL da haste do objeto-modelo, onde são mostradas as forças que agem na haste. 


Aplicando a Primeira Condição de Equilíbrio (Equilíbrio de Translação) e Segunda 
Condição de Equilíbrio (Equilíbrio de Rotação em torno do ponto O da haste), temos: 


F,+F,=F (6.a) (1º Cond. de Equilíbrio) 
F,L, = F,L, (6.b) (2º Cond. de Equilíbrio) 


Sendo, L = (L + L,) o comprimento da haste, isolando L,= L — L., substituindo na 


Eq. (6.b): 
FL, =F;(L—Ly)>(F,+F,)L, =F;L (6.0) 
Logo, Eq. (6.a) em Eq. (6.0): 
FL =F(L-L) > B=(2)F (6d) 
Analogamente, 


F, = (=) F (6.e) 


L> F Ly F F 
K=>T— ,X =T>,X=D— 6.f 
KG LKR nm (6.1) 
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Substituindo Eq. (6.1) na Eq. (5), 


E qe E E A o 
D+ Ps ke) = (6.9) 


Reordenando, 


ke = Dn (7) 


(2) EH) É 


3.2 Contrastação do Modelo Teórico 


Para validar o modelo teórico proposto, contrastaremos os valores de k, que prediz o 
modelo segundo a equação (7) com as evidências experimentais. Para isto, primeiramente, 
são determinadas as constantes elásticas individuais das molas 1 e 2, usando os dados 
coletados de massa da carga (arruelas inseridas num gancho) e da elongação das molas 
medidas com uma balança e uma régua milimétrica, respectivamente. A incerteza da 
medida da massa e da elongação foi de 0,01g e 0,002m, nessa ordem. Posteriormente, 
usando o software origin são plotados os dados da Força F (Peso das “arruelas + gancho” 
= mg) — versus - elongação e via ajuste linear é obtido a constante elástica de cada mola 
através do coeficiente angular, obtendo-se: k = (7,2 + 0,2) N'me k,= (12,8 + 0,2) N/m. 

Posteriormente, as molas 1 e 2 foram suspensas em um suporte horizontal e 
conectadas em seu extremo inferior por uma fina haste de comprimento L (8cm). Em 
seguida, foram coletados dados em duas situações: (1) quando a carga é pendurada no 
centro da haste; e (2) quando a carga é pendurada na haste próxima da mola 1 (L = 3L/8). 
No que segue, apresentamos os resultados obtidos. 


3.2.1 Contraste do modelo teórico para o ponto central da haste L =L/2 


Foram coletados os dados para a carga pendurada no centro da haste, configuração 
experimental representada na Figura 4.b, na qual a haste não permanece horizontal, 
sendo a elongação da mola 1 maior da elongação da mola 2. Em consequência, as forças 
restauradoras das molas F, e F, são diferentes. Os dados experimentais de Força e 
elongação equivalente foram plotados (ver Figura 7), da curva de ajuste linear se determina 
k. =(18,7 + 0,4) N/m ao passo que o valor predito pela Eq.(7), deduzida dentro do modelo 


exp 


teórico, é k. = (18,4 + 0,4) N/m. 
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A molal XX 


1 


o mola? k, 
x mola 142 E 


Força (N) 


k 7,2 0,2)N'm 





k, <12,8+0,2)N'm 
ko (18,7 +0,4)N'm 
k.=18,4+0,4)N'm 
0,00 0,02 0,04 0,06 0,08 
Elongação (im) 


Figura 7: Força em função da elongação para as molas 1, 2 e a mola equivalente. Nesta última, 
as molas 1e 2 são associadas em paralelo, com a carga pendurada no centro da haste. A 
predição do modelo para a constante elástica equivalente está denotado por k... 


Considerando a incerteza do prognóstico do modelo, constata-se que a predição e 
o valor estimado experimentalmente se justapõem, o que pode ser considerado como uma 


evidência da validade do modelo teórico para representar o evento investigado. 


3.2.2 Contraste do Modelo Teórico Para o Ponto da Haste L =3L/8 


Em forma similar que na anterior sub Seção, foram coletados os dados para a carga, 
esta vez pendurada no ponto da haste L = 3L/8 (3cm), configuração experimental similar 
a representada na Figura 4.c. Os resultados são mostrados na Figura 8, onde a constante 
elástica das duas molas associadas em paralelo determinada pelo de ajuste linear da 
curva é KoE (15,6 + 0,4) N/m. Nessa configuração o modelo prediz pela Eq.(6), o valor 
de k.= (15,3 + 0,4) N/m. Portanto, podemos inferir que o resultado experimental corrobora 


empiricamente o modelo teórico contrastado. 
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k 7,2 +0,2)N/m 
k,=(12,8+0,2)N'm 
ko (15,6+0,4)Nim 
k.=(15,3+0,4)Ním 





0,00 0,01 0,02 0,03 0,04 0,05 
Elongação (m) 


Figura 8: Força em função da elongação da mola equivalente, com a carga aplicada em 
L =3L/8. A predição do modelo para a constante elástica equivalente está denotado por k.. 


41 CONSIDERAÇÕES FINAIS 


Os resultados encontrados no presente trabalho são de interesse na discussão 
da construção de modelos científicos, e pela sua simplicidade conceitual e baixo custo 
experimental pode ser desenvolvido, dependendo do enfoque, com alunos do Ensino Médio 
e alunos do Ensino Superior da área de exatas. Um guia de referência para a aplicação 
da modelagem científica em atividades de ensino de Física é exposto com detalhe por 
Heidemann, Araujo e Veit (2016). 

A boa escolha do objeto-modelo, as hipóteses e simplificações deste resultará 
preponderante na obtenção do modelo teórico, além de uma dose de intuição, domínio 
conceitual e uso da matemática. A contrastação empírica do modelo teórico tanto da 
associação de molas em série, quanto, a associação de molas em paralelo foi satisfatória, 
sendo que a hipótese dada pelas equações (2) e (5) são validadas indiretamente através das 
equações (4) e (7), respectivamente. A hipótese (2) e (5) pode ser avaliada conceitualmente 
da forma seguinte: 

No sistema representado pela Figura 5.a, a energia mecânica se conserva. Logo: 


1 1 
= ax + -k2X> — Wa (7.a) 
Onde os dois termos da esquerda da igualdade da Eq. (7.a) são a energia potencial 


elástica das molas 1 e 2, que se opõem à deformação. Por sua vez, o termo à direita, 
é o trabalho realizado pelo campo gravitacional (W.) através da carga pendurada de 
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intensidade F para a associação de molas se deslocar x. Como definido acima F=F. é a 
força variável aplicada à mola equivalente. Logo o trabalho W, é dado por: 


ng = So" Fedxe =, “kexedxs =kexê (rp) 

Ou seja, conceitualmente é demostrado que a soma da energia potencial elástica 
das molas 1 e 2 é igual à energia potencial elástica da mola equivalente. Essa hipótese 
de conservação da energia elástica está embutida na equação (7) que é a predição do 
modelo teórico, o qual foi validado com as evidências experimentais. Destacando, que o 
teste foi feito no modelo teórico e não na Teoria Geral usada “Mecânica Newtoniana” que 


se mantém imune. 
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RESUMO: Os primeiros ímãs que a humanidade 
conheceu foi na antiguidade. Fenômenos 
magnéticos são todas as interações observadas 
entre ímãs, desta forma, o termo magnetismo 
passou a ser empregado para estudar os 
fenômenos e as propriedades que envolvem os 
imãs. Um solenoide, quando percorrido por uma 
corrente elétrica, cria um campo magnético no 
seu interior e exterior, o experimento montado 
para a construção deste trabalho foi chamado 
de Trem Magnético, onde a ideia principal foi 
demonstrar o efeito magnético de um ímã em um 
solenoide, sendo que a montagem é feita com 
materiais de relativo acesso. O experimento foi 
construído e apresentado por alunos da disciplina 
de Laboratório de Eletromagnetismo do curso de 
Licenciatura em Física. 

PALAVRAS - CHAVE: Experimento, 
Eletromagnetismo, Campos Magnéticos. 
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ABSTRACT: The first magnets that humanity 
knew was in antiquity. Magnetic phenomena are 
all the interactions observed between magneis, 
so the term magnetism has come to be used 
to study the phenomena and properties thai 
surround magnets. A solenoid, when traversed 
by an electric current, creates a magnetic field 
inside and outside, the experiment set up for the 
construction of this work was called Magnetic 
Train, where the main idea was to demonstrate 
the magnetic effect of a magnet on a solenoid, 
the assembly is made with materials of relative 
access. The experiment was built and presented 
by students of the Electromagnetism Laboratory 
discipline of the Physics Degree course. 
KEYWORDS: Experiment, Electromagnetism, 
Magnetic Fields. 


11 INTRODUÇÃO 


Os primeiros ímãs que a humanidade 
conheceu foi na antiguidade, na Ásia menor, na 
região da Magnésia, onde foram encontradas 
grandes quantidades de um minério de ferro, 
denominado magnetita, em homenagem a 
região que foi encontrado. A magnetita era um 
material que possuía a propriedade de atrair 
objetos de metal, muito mais tarde os chineses 
usaram as propriedades desse material para 
construir as primeiras bússolas, que possuíam a 
finalidade de orientação nas viagens marítimas 
(MACHADO, 2013). 

Durante diversos 


vários séculos, 


experimentos foram feitos para se obter uma 
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explicação do fenômeno que consistia no alinhamento da agulha de uma bússola ao coloca- 
la próxima de um ímã. Em um experimento realizado por Hans Christian Orsted, em 1820, 
foi verificado que quando uma corrente elétrica passa por um fio que está paralelo a agulha 
de uma bússola próxima a este fio, a agulha sofre um desvio na direção perpendicular ao 
norte-sul. Sendo assim, deduz-se que a corrente elétrica gera um campo magnético que 
interage com o campo magnético da bússola, essa ideia foi relacionada na Lei de Biot- 
Savart (MACHADO, 2013). 

Na Física, denominam-se fenômenos magnéticos todas as interações observadas 
entre ímãs, desta forma, o termo magnetismo é empregado para estudar os fenômenos e 
as propriedades que envolvem os ímãs. Descreve-se abaixo algumas propriedades desses 
materiais: 


| —- Nas primeiras observações realizadas a respeito dos ímãs, foi possível observar 
que eles tinham a capacidade de interagir entre si e também atraíam outros metais. 
Através de experimentos, percebeu-se que esses materiais produziam um campo 
ao seu redor, denominado campo magnético, onde as extremidades desse campo 
foram chamadas de polos do ímã, ou polos magnéticos, designados de polo norte 
e polo sul. 


|| — Colocando um ímã suspenso, percebeu-se que ele sempre se posiciona em 
direção ao norte-sul geográfico. Convencionou-se que polo norte do ímã aponta 
para o polo sul geográfico, e o polo sul do ímã aponta para o norte geográfico. 


||| — Os ímãs exercem entre si forças de ação mútua de atração e repulsão, conforme 
a posição em que são postos um em frente ao outro. Polos iguais se repelem e polos 
diferentes se atraem (SILVA, 2019). 

Uma bobina helicoidal, composta por um fio condutor, formada por espiras circulares 
muito próximas, recebe o nome de solenoide. Um solenoide, quando percorrido por uma 
corrente elétrica, cria um campo magnético no seu interior e exterior, apresentando uma 
configuração de campo magnético semelhante ao de um ímã em barra, então, diz-se que 
ele constitui um eletroímã, ou seja, um ímã obtido por meio de corrente elétrica (HALLIDAY; 
RESNICK; WALKER, 2013). 

O campo magnético do solenoide é a soma vetorial dos campos produzidos 
pelas espiras. No caso de pontos muitos próximos de uma espira, o fio se comporta 
magneticamente como um fio retilíneo, e as linhas do campo são praticamente círculos 
concêntricos. As linhas de campo apresentam pequenos espaçamentos no interior do 
solenoide, o que indica que o campo magnético nessa região é mais intenso. Do lado 
de fora do solenoide, as linhas de campo são mais espaçadas indicando que o campo é 
menos intenso (HALLIDAY; RESNICK; WALKER, 2013). 

A orientação do campo magnético no interior do solenoide é dada pela regra da 
mão direita, segurando o solenoide com a mão direita, os dedos são apontados no sentido 
da corrente, o polegar estendido mostra a orientação do campo magnético (HALLIDAY; 
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RESNICK; WALKER, 2013). 

O experimento montado foi chamado de Trem Magnético simples, onde a ideia 
principal do experimento foi demonstrar o princípio de funcionamento dos trens magnéticos 
reais que são utilizados nas grandes cidades. Em que são utilizados grandes eletroimãs 
para fazer o trem flutuar, deslocar-se e fazer curvas através da força magnética, fazendo 
com que haja diminuição de atrito e um alto desempenho em sua velocidade. Foi idealizado 
e patenteado por Alfred Zehden em 1902 (MONTEIRO, 2018). 

Para a construção deste trabalho, foi realizada uma pesquisa sobre experimentos 
que envolviam fenômenos eletromagnéticos, os quais fossem montados com materiais de 
relativo acesso. O experimento foi construído e apresentado por alunos da disciplina de 
Laboratório de Eletromagnetismo do curso de Licenciatura em Física. 


21 METODOLOGIA 


O experimento realizado para a construção deste trabalho foi baseado em 
experimentos semelhantes encontrados no YouTube. O experimento nomeado de “Trem 
magnético simples” é uma demonstração dos efeitos magnéticos de um ímã em um 
solenoide, sendo que a montagem é feita com materiais de relativo acesso e de fácil 
confecção. O referencial básico para a elaboração do aparato experimental foi um vídeo 
específico do canal Manual do Mundo, o link para o acesso do mesmo está disponível nas 
referências. 

Para a compreensão do funcionamento e dos conceitos fundamentais do fenômeno 
físico gerado pelo experimento, foram feitas pesquisas em livros, artigos, sites e revistas, 
e também, pelas orientações do professor ministrante da disciplina, este projeto foi 
desenvolvido na disciplina de Laboratório de Eletromagnetismo, do curso de Licenciatura 
em Física do Instituto Federal do Piauí (IFPI) — Campus Picos. Ressalta-se que parte da 
montagem foi realizada no laboratório de Física da referida instituição. 

No Quadro 1 são ilustradas, em ordem cronológica, as atividades realizadas para a 


elaboração do projeto, construção e execução do experimento. 
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| Data | Registro de atividades 
01/06/2019 | Pesquisa em sites, artigos e vídeos no YouTube. 


06/06/2019 | Seleção dos experimentos para análise. 
Análise dos experimentos e escolha do tema para o projeto. 


11/06/2019 | Leituras e pesquisas acerca do tema escolhido. 


14/06/2019 | Levantamento de dados e materiais para a confecção do 
experimento. 


19/06/2019 | Confecção do experimento. 
24/07/2019 | Apresentação do projeto. 


Quadro 1 — Cronograma de atividades. 





Fonte: Própria dos autores/2019. 


Para a construção do “Trem magnético simples” foram selecionados os materiais 
que melhor se adequavam com a proposta inicial: a elaboração de um experimento caseiro 
simples. Os componentes utilizados para a confecção do experimento são de fácil obtenção, 
os quais podem ser encontrados em lojas de construção, lojas de eletroeletrônicos e em 
supermercados. Para o funcionamento do trem magnético são necessários os seguintes 


materiais: 


* 1 pilha AAA (diâmetro de 1 mm); 

* | 2imãs (circulares) de neodímio, com diâmetros de 1,5 mm; 

* | 1 rolo de arame de estanho de 1 mm de espessura; 

. 1 cilindro para modelar o estanho, com um diâmetro de aproximadamente 3mm; 


* || 1 tábua de madeira com 40 cm de comprimento e 30 cm de largura. 


A primeira parte do experimento foi a escolha de uma base de suporte para o arame 
de estanho, para isto, foi utilizada uma tábua de madeira (a escolha desse material foi 
devida a sua rigidez). 

Para a segunda parte, tem-se a montagem dos trilhos do trem. Primeiramente é 
necessário modelar o arame de estanho em um formato helicoidal, como um solenoide, 
para isso, usou-se um cilindro que possuía um diâmetro maior do que o diâmetro do ímã, 
cerca de 3 mm, para manter o formato helicoidal, todas as voltas foram colocadas o mais 
próximo possível umas das outras, logo depois de concluído os trilhos, fixou-se na base de 
madeira, como mostrado na figura a seguir (Figura 1). 
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Figura 1 - Montagem do trem magnético. 


Fonte: Própria dos autores/2019. 


As polaridades dos ímãs foram definidas usando uma bússola. Para identificá-los, 
pintou-se o polo norte dos dois ímãs de vermelho, e então foram fixados nas extremidades 
(positivo e negativo) da pilha de forma que os polos de mesmo nome ficassem nas 
extremidades. 

Após a configuração da pilha e dos ímãs, colocou-se o conjunto (pilha e ímãs) dentro 
do solenoide, modelado com o cilindro, e para o imã iniciar o movimento, foi necessário um 
primeiro impulso. O projeto foi apresentado seguindo os passos: 


* | Apresentação do contexto histórico dos trens eletromagnéticos e do experimen- 
to, seguida da abordagem dos conceitos físicos acerca dos mesmos. 


* | Execução do experimento 


* | Discussão sobre a explicação física da atividade experimental com o professor 
e os colegas de turma em sala de aula. 


31 RESULTADOS E DISCUSSÕES 


Após iniciar o movimento do conjunto pilha e imãs (trem) dentro do formato modelado 
de estanho (trilhos), o conjunto pilha e ímãs percorreu todo o comprimento da pista. Este 
movimento se justifica em virtude da interação magnética entre o campo magnético do 
eletrofmã formado pelo solenoide-pilha, e os campos magnéticos dos ímãs. 

De acordo com o posicionamento predefinido dos ímãs na pilha, a força magnética 
proveniente da interação do eletroífmã com os dois ímãs possui a mesma direção e sentido, 
a qual fez mover o conjunto pilha e imãs. A figura 2 ilustra as posições dos ímãs, e a direção 
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e o sentido da força magnética. 





Figura 2 - Posição dos ímãs e sentido da força magnética. 


Fonte: https://encrypted-tbn0O.gstatic.com/ 
images?q=tbn:ANd9GcRPrpq7hcuADL5SWbVwMmik46Fuq94U9P9Q205FGlldVbdy4dd Jw 


Desta forma, as posições dos ímãs nas quais os polos de mesmo nome ficam nas 
extremidades, ou seja, NS e SN causa uma interação de atração, nesse caso do lado 
direito, sendo que no lado esquerdo a interação é de repulsão, fazendo com que a força 
magnética dos dois lados tenha o mesmo sentido, na qual ocorreu a geração do movimento. 

É importante observar que o conjunto pilha e ímãs dentro do solenoide de estanho 
só fecham um circuito elétrico, e consequentemente o eletroímã, porque os imãs funcionam 
como condutores, então a pilha gera uma corrente elétrica que passa pelo imã e em seguida 
no solenoide. Quando o conjunto pilha e ímãs são arrastados, a distribuição espacial da 
corrente no enrolamento acompanha o conjunto, garantindo que a força magnética continue 
a existir sobre ele. 


41 CONCLUSÕES 


Após a apresentação do projeto foi feita uma reflexão analítica sobre a construção 
histórica do experimento, como também sobre os conceitos físicos envolvidos no 
funcionamento e suas interações. Por meio desse projeto, pôde-se ter uma noção 
histórica do surgimento dos conceitos físicos sobre o magnetismo, foi possível observar as 
possíveis interações entre tais conceitos, e principalmente, foi praticável a construção de 
um experimento que simula um caso prático real de um trem magnético, o qual relaciona 
todos esses conceitos. Ao analisarmos o experimento, pode-se identificar as seguintes 
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leis envolvidas: Lei das interações dos polos magnéticos; Lei de Biot-Savart, que são 


fundamentais no estudo do eletromagnetismo. 
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RESUMO: Partindo dos conteúdos ensinados 
no curso de graduação em Física relacionados 
a dinâmica de rotações, o trabalho explica o 
funcionamento do giroscópio, pião, tip-top e 
spinner, os quais, apesar de serem normalmente 
utilizados para recreação, possuem também 
aplicações educacionais e outras utilidades muito 
interessantes de serem exploradas na área de 
Ensino de Física. 

PALAVRAS - CHAVE: Brinquedos. Rotação. 
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ABSTRACT: The first magnets that humanity 
knew was in antiquity. Magnetic phenomena are 
all the interactions observed between magneis, 
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APLICAÇÕES FÍSICAS 


so the term magnetism has come to be used 
to study the phenomena and properties that 
surround magneis. A solenoid, when traversed 
by an electric current, creates a magnetic field 
inside and outside, the experiment set up for the 
construction of this work was called Magnetic 
Train, where the main idea was to demonstrate 
the magnetic effect of a magnet on a solenoid, 
the assembly is made with materials of relative 
access. The experiment was built and presented 
by students of the Electromagnetism Laboratory 
discipline of the Physics Degree course. 
KEYWORDS: Experiment, Electromagnetism, 
Magnetic Fields. 


11 INTRODUÇÃO 


Platão já estudava os movimentos 
celestes, partindo de uma concepção de 
geocentrismo, considerando estes movimentos 
como circulares e harmônicos. Surgiu, assim, 
um grande interesse em estudar objetos que 
rotacionassem ou que apresentavam ciclos. 
Com Newton, no século XVII, toda a Mecânica 
foi revolucionada, uma vez surgiram foram 
cunhados os conceitos de Força, Torque, 
Momento Angular e Linear, utilizando estas 
ferramentas matemáticas que tomariam a 
forma, posteriormente, de notação vetorial. Com 
esta matematização da natureza, foi possível 
estudar com precisão os movimentos e ações 


descritos por objetos em rotação. 
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21 FIDGET SPINNERS 


Os chamados fidget spinners ou hand spinners tiveram origem em 1996 quando 
Catherine Hettinger teve a ideia de criar um brinquedo para crianças com distúrbios de 
ansiedade. Porém, não conseguindo lançar o produto em grande escala, em 2005 ela 
desistiu da patente. No início do ano de 2017 um brinquedo muito similar ao projetado por 
Catherine foi lançado e se espalhou em grande escala por todo o mundo, passando a ser 
fabricado por diversas marcas de brinquedos infantis. 

Basicamente o fidget spinner é um brinquedo em forma de hélice (imagem 1) que 
tem geralmente de duas a três pás (braços). No centro há um rolamento e na extremidade 
de cada braço há um elemento metálico que torna a massa da estrutura do spinner, 
geralmente plástica, desprezível. 

Para utilizar o brinquedo, geralmente segura-se o rolamento utilizando os dedos 
polegar e indicador de uma mão e, com a outra, aplica-se uma força em um dos braços 
do spinner, fazendo-o girar. Com essa rotação surgirá um momento angular no sentido do 
torque e, ao movimentar o spinner, o momento é sentido como uma força resistente ao 


movimento da translação do brinquedo. 





Imagem 1 - Modelo de fidget spinner. 
Fonte: DEPOSITPHOTOS, 2017 (adaptado pelos autores) 


2.1 Funcionamento 


Uma das principais diferenças entre os spinners é a velocidade e o tempo de duração 
da rotação após aplicar o torque. Considerando um spinner (imagem 2) de massa k com n 


“braços” o seu momento de inércia | será: 


Física: O Universo e os Fenômenos Naturais Capítulo 7 HE 


[= Myry + Moy1;?+...+Mary 


onde M é a massa do “braço” e ro seu comprimento. Em um caso onde todos os 
braços são iguais, tem-se que, 
mm ap À É 
[ = nMré., 


onde / é o momento de inércia do spinner. 


Eixo de Rotação 








Imagem 2 - Grandezas presentes na rotação do spinner. 
Fonte: DEPOSITPHOTOS (adaptado) 


Vale ressaltar que quando o spinner está rotacionando de forma equilibrada no 
espaço, o sentido do torque e do momento angular são o do eixo de rotação do spinner, 
sendo que este eixo é o que passa pelo centro de massa e é normal ao plano do movimento 


de rotação do brinquedo. 


31 PIÕES 


A origem do pião é incerta, entretanto existem evidências de que este existe desde 
o ano 4000 a.C. graças a registros arqueológicos na antiga Mesopotâmia. O pião foi trazido 
ao Brasil pelos portugueses e se espalhou rapidamente pelo território nacional, sendo 
muito popular até os dias hodiernos. 

O pião consiste em um brinquedo normalmente feito de madeira, que remete à 


clássica e popular brincadeira de enrolar um cordão em sua ponta e puxar rapidamente 
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fazendo com que este entre em estado de rotação, ficando em pé até o fim do movimento. 
Apesar de seu simples funcionamento, possui uma grande importância para a mecânica 
clássica. 


3.2 Funcionamento 


O pião mantém-se em pé enquanto gira devido à conservação do momento angular 
(L), que tem como direção o eixo de rotação do brinquedo. Para que a direção do momento 
angular mude, é necessário que haja um agente externo que faça mudar sua direção, 
nesse caso, um torque (T). 

Quando o movimento de rotação do pião se inicia (imagem 3), seu eixo e seu 
momento angular têm direção vertical. Conforme continua a girar, este perde velocidade 
devido ao atrito com o chão. Essa queda na velocidade gera outro torque, que faz com que 
ele se incline, criando um certo ângulo 6 entre o chão e o eixo de rotação. Esse movimento 
é chamado de precessão do pião. Neste instante, o momento angular sofre a mesma 


inclinação que o eixo de rotação, fazendo com que ele e o torque se dissipem pelo atrito. 
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Imagem 3 - Esquema do movimento de precessão do pião. 


Fonte: NEVES, M. C. D.; PEREIRA, R. F., 2006, p. 91 


O torque (T) produzido é paralelo ao solo e perpendicular à força peso e ao momento 
angular. É importante ressaltar que ele não vai alterar a intensidade do momento angular, 


apenas a sua direção. O pião continua sob a ação do torque até que se deite no chão (fim 
do movimento). Por definição: 


T=rxm.gousejaT = r.m.g.cosoO 
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Ainda, considerando dt como um pequeno intervalo de tempo e dL um infinitésimo 


de L, tem-se: 


dL 
T=— > dL=Tdt 
dt 


Essa equação descreve um cone (imagem 4), o que explica o fato de a precessão 
do peão ocorrer nessa forma. 


Tdt=dl 


e 4 Tot=i+di 


Imagem 4 - Esquema de ação do torque na precessão do pião. 
Fonte: NEVES, M. C. D.; PEREIRA, R. F., 2006, p. 92 


3.3 Pião Tip-Top 
Um outro tipo de pião bastante conhecido, além do tradicional, é o tip-top (Imagens 
5 e 6). Por ter uma base arredondada ao invés de uma agulha, é mais fácil colocá-lo em 


rotação. Quando posto em rotação, devido ao atrito de sua base com o chão, o tip-top tende 
a se inverter. 





Imagem 5 - Tip-top em funcionamento. 
Fonte: NEVES, M. C. D.; PEREIRA, R. F., 2006, p. 94 
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Ao girar o tip-top, a força de atrito que surge é perpendicular ao plano em que está 
girando, exatamente em seu ponto de apoio. Desse ponto de apoio até seu centro de 
massa há um “braço de alavanca” que causa um torque que permite que o pião se inverta. 
Assim, quanto mais áspera a superfície, mais rápido o pião se inverterá. 

É possível montar uma pião tip-top utilizando 4 bolinhas de gude do mesmo 
tamanho e durepoxi. Basta colar três bolinhas juntas em de triângulo e adicionar a quarta 
no incentro da figura, formando uma pirâmide. Ao colocar o tip-top numa superfície apoiado 
na base triangular, deve-se rotacioná-lo pela bolinha restante. Após certo tempo o tip-top 
se inverterá e começará a rotacionar apoiado nesta. 





Imagem 6 - Tip-top feito com bolinhas de gude e em repouso em rotação. 


Fonte: Elaborada pelo autor 


4 | GIROSCÓPIO 


O nome “giroscópio” (imagem 7) foi criado pelo físico francês Léon Foucault, o 
qual utilizou o princípio de conservação do momento angular para criar um experimento 
relativamente simples capaz de provar a existência do movimento de rotação da Terra. 

Um giroscópio é essencialmente formado por um rotor perfeitamente balanceado que 
gira livremente entre seus eixos geométricos perpendiculares entre si que se interceptam 
no centro de massa do rotor. Quando posicionado numa suspensão é denominado de tipo 
Cardan. O giroscópio, quando em rotação, pode adquirir qualquer orientação, porém seu 
centro de massa deve estar sempre fixo no espaço. 
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Gyroscope Spin axis 
frame 






Gimbal Rotor 


Imagem 7 - Giroscópio em suspensão Radan. 
Fonte: WIKIMEDIA, 2017. 


Devido à ação da inércia giroscópica, o rotor tende a conservar sua direção no 
espaço, independentemente dos movimentos aplicados no mesmo. Ou seja, a inércia 
giroscópica é uma propriedade que o giroscópio tem em manter seu centro de massa na 
mesma posição. 

Para entendermos a experiência de Foucault imaginemos um giroscópio na posição 
A (imagem 8) da figura, com eixo do volante apontando para a horizontal, e cuja direção 
deve permanecer constante durante todo o movimento. Se a Terra não girasse, este deveria 
permanecer sempre na horizontal, porém, o observado é que, após algumas horas, quando 
o objeto estiver em B, a direção do eixo agora é a vertical, comprovando a rotação de nosso 
planeta. 





Imagem 8 - Experimento esquematizado de Foucalt. 


Fonte: NUSSENZVEIG, 2002, p.276 
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O experimento original ocorreu em Paris no ano de 1851. Foucalt posicionou um 
giroscópio conhecido como Pêndulo de Foucalt no Panthéon, cuja direção do eixo de 
rotação do volante variava de acordo com esperado para o movimento de rotação da Terra. 

O fato de a direção do eixo de rotação do giroscópio se conservar permite utilizá-lo 
em diversas áreas. Pode ser usado em navegação marítima, aérea e espacial para evitar 
que o veículo se distancie de uma direção pré-estabelecida, tarefa que este efetua com 
grande precisão, uma vez que é minimamente afetado por perturbações eletromagnéticas 
ou pela inclinação lateralNongitudinal do transporte. Possui também aplicações militares, 
sendo utilizado para guiar mísseis teleguiados, além de seu emprego na orientação da 
perfuração de poços de petróleo. 


5 | CONSIDERAÇÕES FINAIS 


O estudo analisou o funcionamento de brinquedos envolvidos com a rotação com 
foco nos princípios físicos (velocidade e aceleração angular, centro de massa, momento de 
inércia, torque, momento angular e sua conservação) que explicam tal movimento. 

O caráter experimental da Física torna necessárias demonstrações experimentais 
para facilitar o aprendizado, tarefa que não é possível ser realizada para todas as áreas 
dessa ciência, uma vez que, na maioria dos casos, faltam recursos para sua consecução. 
Além disso, os estudantes passam uma grande parte do tempo de estudo resolvendo 
problemas distantes de seu cotidiano. A utilização dos brinquedos em sala de aula para 
demonstrar os conceitos físicos previamente citados resolve parcialmente estes problemas, 
aproximando o dia-a-dia dos alunos do estudo da Física, tornando a experiência mais 
prazerosa e gerando resultados que podem consolidar o aprendizado. 
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RESUMO: O Ppresente trabalho tem por 
finalidade o estudo de trajetórias bidimensionais 
de forma integrada, aliando a confecção de uma 
catapulta caseira, em escala, no estilo Trebuchei 
a realização de um estudo aprofundando dos 
conceitos físicos necessários para descrição 
de seu funcionamento, bem como, a dedução 
de funções parabólicas capazes de prever com 
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considerável exatidão a real trajetória de um 
projetil sob efeito, unicamente, da força resistiva 
do ar. A extração de dados experimentais e 
simulações das trajetórias, foram obtidas por 
intermédio do software livre Tracker, ferramenta 
de vídeo-análise voltada para o ensino, e 
linguagem de programação em Fortran 90. 
PALAVRAS - CHAVE: Lançamento de projétil; 
Constante da força de retardo; Trebuchet; 
Tracker, Fortran 90. 


TREBUCHET: EFFECT OF THE AIR 
RESISTANCE ON THE PROJECTILE 
MOVEMENT 


ABSTRACT: The present work has the purpose 
the study of two-dimensional trajectories in an 
integrated manner, combining the building a 
homemade catapult in scale in the Trebuchet 
style and the realization of a deep study of the 
physical concepts necessary for the description 
of its operation, as well, the deduction of parabolic 
functions able to predict with considerable 
accuracy the real trajectory of a projectile under 
the effect of the resistive force of the air, only. The 
extraction of experimental data and trajectory 
simulations, were obtained through the free 
software Tracker, a video-analysis tool aimed 
at teaching, and the Fortran 90 programming 
language. 

KEYWORDS: Projectile launch; Drag coefficient; 
Trebuchet; Tracker, Fortran 90. 
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11 INTRODUÇÃO 


O processo de aprendizagem dos conceitos de Física, pode ser facilitado com a 
experimentação e a produção manual, por parte dos alunos, de modelos que simulem os 
conteúdos propostos em aula. (BORGES, 2006) Neste sentido foi proposto aos alunos do 
Curso de Bacharelado de Física da Universidade Estadual de Ponta Grossa um estudo 
detalhado do deslocamento balístico de um projétil, lançado através de uma catapulta 
clássica denominada Trebuchet. 

O advento do Trebuchet remonta aos povos chineses que, entre os séculos V e III 
AC, desenvolveram um modelo primitivo que dependia exclusivamente da força humana 
no acionamento da arma de cerco, via conjunto de cordas e polias (CHEVEDDEN, 1995; 
SAIMRE, 2006). A partir do século XII DC, especialmente em países islâmicos, bizantinos 
e em áreas mediterrâneas da Europa, há relatos do emprego do Trebuchet operado por 
ação de contra peso no arremesso de seus devastadores projéteis. A França figura como o 
primeiro país europeu a empregar o modelo de acionamento por contra peso (ANNÓSCIA, 
2018). 

Devido a sua simplicidade de construção e emprego, o Trebuchet resistiu bravamente 
frente ao surgimento dos primeiros canhões a base de pólvora e, somente por volta do 
século XV DC, teve seu emprego gradativamente descontinuado devido a novos canhões 
capazes de rivalizar com o poder de destruição da arma de cerco (SAIMRE, 2006). 

Neste contexto histórico, o conhecimento da mecânica teve que ser desenvolvido 
portanto, é justo dizer que o trebuchet contribuiu para o progresso da ciência à sua maneira. 

Aliando os conhecimentos históricos aos princípios básicos de mecânica, busca-se 
a confecção de um modelo caseiro em escala reduzida do Trebuchet. Neste trabalho faz- 
se um estudo didático do funcionamento deste aparato com auxílio do software Tracker 
(PARREIRAL 2018) e linguagem Fortran 90, bem como a comparação das trajetórias 
parabólicas, modeladas para condições ideais (sem forças dissipativas) ou adaptadas, 
unicamente, para presença de um fator da intensidade da força de amortecimento no 
sistema. 

A construção da catapulta pode ser feita de forma caseira e com material reciclado. 
Deve servir de estímulo aos alunos do ensino médio e graduação dos cursos de exatas e 
engenharias a entender a física por trás do funcionamento de uma antiga arma de guerra. 
1.1 Trebuchet 


O Trebuchet consiste em um conjunto composto por duas partes básicas: a estrutura 
de sustentação e as partes móveis. O funcionamento resume-se ao movimento das partes 
móveis sobre sua estrutura base, permitindo que a alavanca juntamente com o contrapeso 
arremesse o projétil através do estilingue. A imagem da figura 1, a seguir, ilustra O arranjo 


descrito. 
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Figura 1 — Arranjo esquematizado de um Trebucheti 





Conforme descrito acima, o projeto em miniatura baseia-se no modelo histórico, 
apresentando algumas alterações necessárias devido às limitações nos materiais 
disponíveis. Podemos visualizar as dimensões do protótipo na figura 2 e ele confeccionado 
na figura 3, de acordo com as especificações anteriormente mencionadas. 





Figura 2 — Esquema das dimensões do protótipo. 
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Figura 3 — Imagem contendo o modelo em escala do Trebuchet construído de acordo com as 
especificações propostos. 


O projétil utilizado nos lançamentos foi um corpo de prova com dimensões 
comparáveis a um geoide, sendo confeccionado artesanalmente com massa adesiva epóxi 
e com aproximadamente 50 q. 

A obtenção dos valores correspondentes à velocidade inicial (v.), tempo de voo 


(t), alcance máximo (A |), ângulo de arremesso (e), distância do projétil ao centro da 


máx 
alavanca (1), entre outros parâmetros; foram aferidos através da recordação em vídeo para 
cada lançamento e sua posterior análise por meio do software livre Tracker. A região para 
lançamentos foi limitada para um alcance máximo de 11,5 m onde a posição inicial (x,) 
corresponde à média das posições no instante em que ocorre o desprendimento entre o 
corpo de prova e o aparato. A construção gráfica e análise das trajetórias — nas condições 
sem arrasto e com arrasto — foram intermediadas por um algoritmo escrito em linguagem 


Fortran 90. 


1.2 Trebuchet e o movimento de projétil 


A energia envolvida no processo é simplesmente a energia potencial gravitacional U 
; dada pela equação (1) a seguir: 


U=mgh (1) 


com uma aceleração gravitacional g e de módulo constante de aproximadamente 
9,79 m/sº — Ponta Grossa, Paraná. 
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Fica evidente que ao acionarmos o Trebuchet, no instante em que o contra peso 
é liberado, um movimento semelhante ao de queda livre se iniciará, convertendo toda U 
acumulada no sistema em energia cinética de rotação, cinética de translação e em energias 
dissipativas — atrito, calor, som, etc., resultando no arremesso do projétil que se encontra 
na outra extremidade do artefato (HALLIDAY, 2009) 

De todas as peças que compõem o Trebuchet, pode-se dizer que a alavanca é 
aquela que desempenha um papel chave em seu funcionamento. O torque resultante Tres 
sobre o sistema (HALLIDAY, 2009; NUSSENZVEIG, 2002), segue a equação abaixo: 


Tres = Fprsenb (2) 


Conforme o conjunto realiza um movimento de rotação devido ao torque resultante, 
o corpo de prova se desloca de sua posição inicial até atingir o ponto de sua liberação do 
estilingue. Durante o deslocamento, há uma ampliação na distância r que separa o projétil 
do eixo de rotação bem como um aumento gradativo no momento de inercia total 1.,, € 
no módulo da velocidade angular w na extremidade composta pela alavanca, estilingue e 
corpo de prova. Como resultado, parte da energia potencial gravitacional U é convertida em 
energia cinética de rotação K. Aimagem da figura 4, nos mostra a realização do movimento, 


de forma que K acaba sendo escrita como: 


1 
K, = 3 Irotaiw? (3) 


Onde temos o momento de inercial total Urotai) é módulo da velocidade angular (w). 
O momento de inércia (NUSSENZVEIG, 2002) do conjunto no instante do lançamento será 


descrito por: 


Irotal e Iuaste F Epotai 


1 2 2 
Irotal = 12 MT3 + mr3 (4) 


com M à massa da haste e m à massa do corpo de prova. 
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Figura 4 — Representação das posições do projétil durante o lançamento e a variação da distância até 
seu eixo de rotação. 


Após a rotação da alavanca temos a liberação do objeto pelo estilingue, 
arremessando-o com uma velocidade inicial v, e com uma trajetória parabólica. A razão por 
trás de tal funcionamento reside no fato de que, parte da energia potencial gravitacional é 
transformada — além da energia cinética de rotação — em cinética de translação, Kr = mos, 
apresentando valor máximo no instante t, da liberação do corpo (NUSSENZVEIG, 2002). 

Numa primeira análise, se executássemos o arremesso sob condições ideais, sem 
atrito e somente a gravidade agindo sobre o corpo, obteriamos uma função parabólica 
(MARION, 1988) de y em função de x que pode ser descrita por: 


9 


W)=Y+xtan0-x—-—— — 
YO) = 2vê cos? 8 


(5) 

Um olhar atento nos revela que a equação (5) corresponde a uma função do 2º grau 
onde o coeficiente a corresponde aos valores 9/2v; cos? 6, b equivale à tan O e c associa- 
se a Y,. Embora válida, a mesma desconsidera a resistência oferecida pelo ar — ou arrasto 
— durante o deslocamento, responsável por reduzir os módulos de (x) e vo. 

Se desejamos modelar o comportamento de um corpo em trajetória parabólica sob 
efeito de uma força resistiva e, considerando somente a influência do arrasto (CHANG, 
2019), podemos assumir que o módulo da força resistiva FrEs será proporcional ao módulo 


da velocidade v do projetil conforme abaixo: 
Fres « v” (6) 


Quando objetos pequenos se deslocando com velocidades que sejam 
aproximadamente menores que 24 m/s, o fator »n assume um valor de n = 1 (MARION, 
1988). Logo, o módulo da força resistiva que atua sobre o corpo de prova será: 


FrEs = —mkv (7) 
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Consequentemente, devemos reescrever a equação (5) para que a constante positiva 
k possa ser levada em consideração, e que mais adiante, possibilitará uma modelagem 
com exatidão do movimento de um corpo sujeito ao arrasto, após seu arremesso. Assim 
sendo, consideremos o esboço do movimento de projétil sob efeito de uma força resistiva 
através da figura 5. 


y V, V 
e 0 V. Aceleração Gravitacional 


g 


Força Resistiva 


Fres =mk V 
Yo 


X 


Figura 5 — Representação, em coordenadas retangulares, da trajetória parabólica de um corpo 
influenciado pela força de arrasto. 


A expressão que satisfaz o modelo de Trebuchet analisado é (MARION, 1988): 


+ (ecoça) (o) (8) 


A equação (8) consiste na solução geral para trajetórias parabólicas de projéteis 


kx 


Vo cos O 











Y(x) = Yo +Elnl1 — 


sob efeito de uma força resistiva do tipo Fres = —-mkv. Os parâmetros expressos em (8) 
podem ser identificados como sendo: equação parabólica da posição vertical em função 
da posição horizontal do projétil (y[x]), altura inicial do projétil (9), módulo da aceleração 
da gravidade local (k), constante associada à intensidade da força de amortecimento (k), 
módulo da velocidade inicial do projétil (v,), posição horizontal do projétil (x) e o ângulo de 


arremesso do projétil (e). 


1.3 Vídeo-análise do movimento e linguagem fortran 90 


A obtenção de parâmetros físicos para o estudo das propriedades mecânicas do 
Trebuchet ocorreu por intermédio do software de vídeo-análise Tracker (PARREIRAL, 2018) 
e a linguagem de programação Fortran 90. (ANAGNOSTOPOULOS, 2016) O primeiro é 
capaz de fornecer dados experimentais atrelados às trajetórias de corpos em movimento 
enquanto que, o segundo, permite o desenvolvimento de programas executáveis capazes 
de compilar e manipular dados de acordo com o que se deseja analisar. 

O software Tracker apresenta uma barra de ferramenta poderosa contendo 
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comandos de ajuste e análise, a saber: eixo de coordenadas — permite orientar o sistema 
através de coordenadas retangulares definido em unidades de pixel, além de comportar o 
ajuste do ângulo formado entre o eixo e o plano do vídeo —, fita de calibração — possibilita 
que o software calibre o tamanho dos objetos em vídeo com suas dimensões reais -, 
bastão de medição — permite a extração precisa das dimensões dos objetos contidos em 
vídeo, estando dependente da utilização correta da fita de calibração —, ponto de massa — 
possibilita a inclusão da massa e rastreamento da trajetória do objeto retornando valores de 
grandezas físicas associadas à posição e massa -, entre outras opções. 

A imagem da figura 6 exemplifica a funcionalidade das ferramentas eixo de 
coordenadas, fita de calibração, bastão de medição, ponto de massa, visualização de 
trajetória; disponibilizadas pelo software Tracker (PARREIRAL, 2018). 





Figura 6 — Imagem extraída da análise do Teste nº5 do corpo de prova. 


Por intermédio da linguagem em Fortran 90 (ANAGNOSTOPOULOS, 2016) 
desenvolveu-se um algoritmo contento os parâmetros experimentais específicos do corpo 
de prova. Desta forma, obteve-se uma comparação gráfica entre as trajetórias descritas 
nas condições sem atrito (ideal) e na presença, somente, da ação dissipativa da resistência 
do ar (intermediada pelo fator da força de amortecimento). 

A tabela abaixo apresenta os comandos necessários para o desenvolvimento do 
algoritmo já mencionado. 
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program Projetil 


implicit none 


real :: y 

write(*,)y() 

end program Projetil 
real function y() 


implicit none 


character(20),parameter :: inputfile ="Dados A” 
character(20),parameter :: outputfile ="Resultados A” 
integer, parameter :: n=104 


integer :: 


real :: x(n),y1,y2 
real,parameter :: A1=0.151,B1=1.41,C1=1.362 
real,parameter :: k=0.175128387,hn=1.362,9=9.79,V0=9.879,A=54.84 


open(1 ,file=inputfile) 


open(2,file=outputfile) 


do i=1,n 

read(1,*)x(i) 

enddo 

do i=1,n 

y1=C1+B1*(x(1))-A1*(x(1)**2) 
y2=h+((g/k**2)*log(1-x(i)*k/VO*cos(A)))+((x(i)/VO*cos(A)*k“*2)*(VO*sin(A))) 
+((x(i)/VO*cos(A)*k**2)*(g/k)) 

write(2,*)x(i),y 1,y2 


enddo 


close(2) 


end function 





Tabela 1 — Programa Corpo de Prova em linguagem Fortran 90. A função denota a trajetória 
sem atrito enquanto que, a função y2, descreve a ação resistiva do ar. 
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2 | RESULTADOS E DISCUSSÕES 


Podemos condensar na tabela que segue os valores médios provenientes dos 
respectivos testes executados com o corpo de prova. Os testes resumiram-se a cinco 


lançamentos do projétil, procedimento mencionado anteriormente. 


ALCANCE MÁXIMO Ay (M) 8,558 
ÂNGULO DE LANÇAMENTO O (graus) 
ACELERAÇÃO GRAVITACIONAL g (m/s?) 





Tabela 2 — Valores médios dos parâmetros experimentais extraídos com auxílio da ferramenta 
Tracker. 


Com auxílio das informações contidas na tabela acima e da equação 5, podemos 
construir uma função parabólica da suposta trajetória do corpo de prova na ausência da 
ação resistiva do ar. Portanto, considere abaixo: 


Ysa = —(0,151)x? + (1,420)x + (1,362) (9) 


Onde a notação y.., representa à função matemática do corpo de prova sem influência 
do atrito. Alternativamente, em posse da equação 8 que prevê a influência do arrasto sobre 
o corpo de prova — denotado por y,., — e, limitando-a no intervalo O < x < 8,558, realizou- 
se uma manipulação numérica do módulo da constante k, de maneira que a função se 
aproximasse de zero (Yca = 0) no instante X = 8,558, valor correspondente ao alcance 
máximo de A. Sendo assim, a função y,, utilizada na manipulação matemática já descrita, 
encontrasse explícita abaixo: 
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-1,362+(25)1 1- aid | 
Yca = À, (Kka)?) “|” (9,879)co s(54,84) o 
ta n(54,84) A (9,79) 10) 
” ka (9, 879)co s(54, 84) 


Conforme o raciocínio prévio, o ajuste para constante K capaz de descrever 
satisfatoriamente a trajetória experimental do projétil equivale a K = 0,175128387 — ou, 
aproximadamente 0,175 —, em unidades de s!. Consequentemente, a comparação entre a 


trajetória sem arrasto (cor azul) e com arrasto (cor vermelha) pode ser observada abaixo: 


— Trajetória sem 
Corpo de Prova Atrito 


— Trajetória com 
Atrito 





Figura 7 — Comparação entre as funções trajetória do corpo de prova. 


Aplicando na equação 7 os parâmetros m = 0,05 kg, vo = 9,879 m/se k = 0,175 5” 
teremos uma força resistiva diretamente proporcional a velocidade do corpo de prova com 
intensidade aproximada de: 


Fes = — 0,086 N (11) 


As trajetórias balísticas contidas na Figura 7 evidenciam uma clara influência do 
arrasto do ar em sistemas físico-mecânicos. Nota-se nesta figura, que para condições ideais 
descritas pela equação (9) o projetil teria um alcance total de aproximadamente 10,264 m. 
Não obstante, ao contabilizarmos os efeitos do arrasto na trajetória efetuada, podemos 
extrair duas informações expressivas: a redução significativa no módulo do alcance total, 
equivalente a 8,558 m, e a elevada acurácia das equações em fornecer o valor do alcance 


máximo teórico (8,5582 m) em relação ao valor experimental (8,558 m). 
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31 CONSIDERAÇÕES FINAIS 


Observa-se que o modelo experimental criado neste trabalho está em excelente 
concordância com a simulação computacional. De sorte que o modelo proposto é adequado 
para o estudo balístico de projéteis. 

Experimentalmente a catapulta, Trebuchet, aqui utilizada é uma ferramenta 
adequado para o estudo e determinação da resistência oferecida pelo ar na propagação 
de um projétil. 
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RESUMO: A experiência a ser relatada se deu 
com duas alunas do Curso Licenciatura em 
Física, que diante da problemática, cujo foco 
consiste no ensino de Física para os alunos 
com surdez do 1º ano do Ensino Médio dos 
Cursos Integrados em Informática e Química 
do Instituto Federal de Educação, Ciência e 
Tecnologia do Amazonas (IFAM), do Campus 
Manaus Centro (CMC). É importante destacar 
que os esses estudantes envolvidos nessa 
pesquisa, frequentavam aulas de reforço da 
disciplina Física, oferecidas pelo Núcleo de 
Atendimento à Pessoas com Necessidades 
Educacionais Específicas (NAPNE), contudo 
para o desenvolvimento dessas aulas, foram 
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AMAZONAS-IFAM/CMC 


elaborados os resumos adaptados e os jogos 
educativos adaptados para auxiliar no processo 
de Aprendizagem. 

PALAVRAS - CHAVE: Ensino de Física; Inclusão; 
Jogos educativos adaptados. 


ABSTRACT: The experience to be reported 
took place with two students from the Physics 
Degree Course, who faced the problem, whose 
focus is on teaching Physics to students with 
deafness in the ist year of High School of the 
Integrated Courses in Informatics and Chemistry 
at the Federal Institute of Amazonas Education, 
Science and Technology (IFAM), Campus 
Manaus Centro (CMC). It is important to highlight 
that these students involved in this research, 
attended classes to reinforce the Physical 
discipline, offered by the Center for Assistance 
to People with Specific Educational Needs 
(NAPNE), however for the development of these 
classes, adapted summaries and educational 
games were elaborated. adapted to assist in the 
Learning process. 

KEYWORDS: Physics teaching; 
Adapted educational games. 


Inclusion; 


11 INTRODUÇÃO E JUSTIFICATIVA 


A experiência a ser relatada se deu 
com duas alunas do Curso Licenciatura em 
Física, que diante da problemática, cujo foco 
consiste no ensino de Física para os alunos 
com surdez do 1º ano do Ensino Médio dos 
Cursos Integrados em Informática e Química 
do Instituto Federal de Educação, Ciência e 
Tecnologia do Amazonas (IFAM), do Campus 


Capítulo 9 


Manaus Centro (CMC). É importante destacar que os esses estudantes envolvidos nessa 
pesquisa, frequentavam aulas de reforço da disciplina Física, oferecidas pelo Núcleo de 
Atendimento à Pessoas com Necessidades Educacionais Específicas (NAPNE), contudo 
para o desenvolvimento dessas aulas, foram elaborados os resumos adaptados e os jogos 
educativos adaptados para auxiliar no processo de Aprendizagem. 

O IFAM ao longo dos seus 110 anos, vem atendendo à demanda de alunos que 
ingressam na instituição procurando se adequar às novas leis de inclusão para a educação. 
Para isso, criou o NAPNE pela Portaria nº 180- GDG/ CEFET - AM/2002, espaço que 
visa promover acessibilidade na instituição com o objetivo de receber estudantes e 
providenciando a adaptação do currículo, conforme a necessidade de cada um. O 
referido Núcleo, disponibiliza atualmente para os alunos com necessidades específicas 
educacionais a acessibilidade arquitetônica, a acessibilidade pedagógica, a acessibilidade 
comunicacional e a acessibilidade atitudinal, além de disponibilizar tradutores/intérpretes 
da Língua Brasileira de Sinais (Libras) para acompanhar estes discentes durantes as suas 
aulas. 

Um ponto importante é que nem sempre o intérprete dispõe de conhecimentos a 
respeito dos conceitos e teorias que os professores ensinam em sala de aula, dessa forma, 
isso acaba sendo um problema, pois estes devem repassar aos alunos o que o professor 
fala em sala de aula, tornando portanto, a tradução para a Libras de forma equivocada. 
Além disso, não existem terminologias em Libras para descrever alguns fenômenos e 
conceitos específicos da Física. Com isso, os intérpretes enfrentam dificuldades com as 
especificidades, já que nem sempre o conhecimento deles, condiz com o conceito que está 
sendo trabalhado, de modo que a tradução para Libras pode ocorrer de maneira distorcida 
(SILVEIRA; SOUSA, 2011). 

A Física é uma disciplina que é considerada complexa, até para pessoas ouvintes, 
sendo esta uma das ciências que lida basicamente com interpretações de fenômenos 
naturais do cotidiano, esses fenômenos são moldados através da matemática, e para 
abordar alguns temas, necessitamos da audição para conhecer os sons graves ou agudos, 
como a Acústica por exemplo. 

Não podemos deixar de destacar que existem conceitos da Física que ainda não 
foram traduzidos para a Libras, devido a este fator torna-se difícil ensinar vários conteúdos 
da Física |, desta forma o professor ou intérprete que acompanha estes discentes, deve 
procurar outros recursos didáticos como os jogos educativos adaptados para facilitar a 
aprendizagem, para ensinar vários conceitos e descrever fórmulas da Física que não são 
compreendidos de forma satisfatória por estes discentes, visto que a forma de ensino 
tradicional não contribui para a aprendizagem dos mesmos e também porque estes 
possuem uma limitação à falta de audição, que embora não interfira no seu aprendizado, 
mas limita a aquisição de conhecimentos transmitidos oralmente, conforme afirma Zenari 


(2003) “a falta de audição não afeta as capacidades intelectuais, mas limita a possibilidade 
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de aquisição de conhecimentos transmitidos oralmente, prejudicando o desenvolvimento 
do raciocínio abstrato, já que haverá dificuldades em formar conceitos simbólicos, que não 
necessitem da exploração concreta dos objetos”. 

Com base nas considerações e problemas apontados no ensino de Física para 
discentes surdos, propõe-se utilizar uma metodologia de ensino diferenciada, com o 
uso dos jogos educativos adaptados, que foram confeccionados durante o projeto, com 
o objetivo de facilitar a compreensão dos assuntos da Física |, estimular a curiosidade 
desses discentes surdos e fazer um acompanhamento do rendimento dos mesmos nesse 
componente curricular. Nosso público-alvo, foram discentes surdos do IFAM-CMC, que 
durante a execução deste projeto contou coma participação de cinco discentes surdos 
matriculados no Ensino Médio Integrado sendo: cinco alunos surdos matriculados no 
primeiro ano, mas em diferentes cursos, três discentes surdos matriculados no curso 
Integrado em Informática, dois discentes surdos matriculados no curso Integrado em 


Química. 


21 OBJETIVOS 


2.1 Objetivo Geral 


Propor uma metodologia de ensino diferenciada, onde através da utilização de 
jogos educativos adaptados, possamos melhorar o ensino de Física | para estudantes com 
surdez do IFAM-CMC. 


2.2 Objetivos Específicos 


* | Com base em um levantamento dos conceitos da Física que não estão contex- 
tualizados e que são ensinados de forma muito teórica, a partir desses temas, 
elaborar jogos educativos para facilitar a aprendizagem dos discentes; 


* Desenvolver uma metodologia de ensino diferenciada, associada a utilização 
de recursos “jogos educativos adaptados” para o ensino de temas da Física 
| disponibilizando uma linguagem simplificada, imagens, termos em Libras e 
utilização de palavras-chaves para facilitar a compreensão e aprendizagem dos 
discentes. 


31 MÉTODO 


Primeiramente foi feita uma pesquisa, onde o objetivo era saber quais os assuntos 
do componente curricular Física 1 que estes discentes mais apresentavam dificuldades. 
Posteriormente, foram elaborados os resumos adaptados sobre vários temas, tais como 
cinemática, movimento em uma dimensão, movimento em duas dimensões e as Leis de 


Newton. Foi levado em consideração que os discentes surdos, compreendem melhor 
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através do aspecto visual, dessa forma, procurou-se elaborar um material com a menor 
quantidade de textos possível e maior quantidade de elementos visuais, como por exemplo 
imagens, diagramas, esquemas e fluxogramas. 

Para a confecção dos jogos educativos adaptados para o ensino dos temas de 
Física, elaboramos resumos de cada conteúdo do bimestre para facilitar a aprendizagem 
dos discentes surdos, todavia, nesses resumos, foram utilizados muitas imagens, exemplos 
do cotidiano, detalhamos às fórmulas de cada conteúdo e ao final de cada resumo, 
propomos um jogo educativo adaptado, contendo a cruzadinha, o caça-palavras ou o jogo 
da memória de acordo com o assunto abordado no resumo. Logo, considerando o conteúdo 
que foi trabalhado nas aulas do reforço de Física na sala de atendimento, associado ao uso 
desses materiais como recursos, nossa ideia é verificar o entendimento dos estudantes nos 
conceitos estudados e as respectivas contribuições na aprendizagem dos mesmos. 

Os exemplos e exercícios que foram utilizados nos resumos adaptados, foram em 
sua maioria tirados da Apostila de Física Anchieta, visto que a mesma apresenta uma 
linguagem simplificada e possui várias imagens. Agora, em relação aos assuntos que não 
eram explicados com uma abordagem visual através dessa apostila, procurou-se pesquisar 
em outras fontes, como livros didáticos e livros digitais. Para confecção deste material 
didático em forma de resumo, utilizamos além das imagens, fluxogramas, que facilitavam 
a compreensão do aluno e permitia a visualização do passo-a-passo da solução de um 
exemplo ou exercício resolvido. Também utilizamos, as palavras-chaves, que foram usadas 
como palavras principais, fazendo, portanto, a associação com o conteúdo abordado. 

Mesmo se considerarmos um aluno comum, este necessita de apoio e 
acompanhamento em suas atividades, pois surgem dúvidas no processo de aprendizagem, 
o que é perfeitamente normal. Em virtude disso, oferecemos apoio aos alunos surdos, 
atendidos pelo projeto para esclarecer os conteúdos e sanar possíveis dúvidas referentes 
aos resumos adaptados e aos jogos educativos adaptados. 

Com base no acompanhamento dos discentes surdos, através das aulas de reforço 
de Física, percebemos a importância de elaborar recursos didáticos adaptados para esse 
público. Para SOUZA (2007, p.111), “Recurso didático é todo material utilizado como auxílio 
no ensino-aprendizagem do conteúdo proposto para ser aplicado, pelo professor, diante de 
seus alunos”. 

Atualmente existe uma variedade de recursos didáticos que são utilizados para 
trabalhar com discentes surdos do Ensino Fundamental, porém existem poucos recursos 
didáticos para se trabalhar com este público no Ensino Médio, público que possui uma 
forma diferente de aprender. 

De posse dessas informações sobre as particularidades dos discentes foram 
elaborados e confeccionados os jogos educativos adaptados sobre os assuntos de Física 
1, como as cruzadinhas, o jogo de caça-palavras e o jogo da memória, sobre temas da 
Física contextualizados com as situações cotidianas destes discentes surdos. É importante 
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destacar que os jogos educativos adaptados foram elaborados e confeccionados de acordo 
com os assuntos do componente curricular Física |, que tinham sido abordados nas aulas 
de reforço, usando a sala de atendimento, que é o local onde são ministradas as aulas 
desse reforço. De acordo com o desempenho dos discentes nas atividades e nos jogos 
educativos adaptados, observou-se que através da inserção desse recurso didático, houve 
uma evolução na aprendizagem desses discentes, pois os mesmos já tinham estudado 
sobre os assuntos da Física | que envolviam os jogos durante as aulas de reforço e sala de 
aula regular, além de ter aumentado a interação com os colegas surdos, ainda estimulou o 


espírito de competição durante os jogos entre os mesmos. 


41 RESULTADO E DISCUSSÃO 


Produzir um recurso didático para o ensino de discentes com surdez não é uma 
tarefa fácil, pois requer muito esforço, trabalho e dedicação por partes dos envolvidos. Se 
tal material for direcionado para um público específico, a tarefa é ainda maior, carecendo 
de um estudo aprofundado sobre o objeto. 

Para a produção dos materiais didáticos de Física para alunos surdos, tivemos que 
seguir alguns passos no desenvolvimento das atividades, que foram necessárias para o 
andamento do trabalho. 


Passo 1: Pesquisa sobre os materiais didáticos produzidos sobre o Ensino de 
Física para discentes surdos 

Inicialmente, foi necessário realizar pesquisas sobre os materiais que tinham sido 
produzidos na mesma linha de pesquisa do projeto, porém ressaltamos que poucos materiais 
foram encontrados, mas, no entanto, alguns artigos sobre trabalhos com alunos deficientes 
auditivos serviram como embasamento teórico, de como proceder nas atividades. 

Lembramos que o trabalho de pesquisa sobre os temas que seriam trabalhados 
nesse projeto, foram realizados no primeiro mês. Também nesse período, foi feito o 
acompanhamento dos discentes surdos nas aulas de reforço de Física, com o objetivo 
de identificar problemas relacionados à forma convencional de ensino e as principais 
dificuldades dos discentes surdos atendidos pelo projeto, conforme é mostrado na tabela 
(1) abaixo: 
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Alunos Curso / ano 
Aluno A INF /1º ano 
Aluno B INF /1º ano 
Aluno C INF /1º ano 
Aluna D QUI /1º ano 
Aluna E QUI /1º ano 





Tabela 1: Discentes participantes do projeto 


Fonte: Acervo pessoal 


Passo 2: Elaboração dos Resumos Adaptados sobre temas da Física 


Para dar início a elaboração dos resumos adaptados que antecediam aos jogos 
educativos adaptados, foi feita uma pesquisa para encontrar um livro didático ou apostila que 
tivesse uma linguagem simplificada, que pudesse auxiliar na elaboração destes resumos 
adaptados. Utilizamos então, a apostila construindo a Física material pertencente ao 
colégio Anchieta, onde é importante destacar que esse material apresenta uma linguagem 
simples e possui várias imagens, desta forma facilitou o nosso trabalho. 

Posteriormente, foram elaborados os resumos adaptados sobre vários temas da 
Física | como cinemática, movimento em uma dimensão, movimento em duas dimensões e 
as Leis de Newton. Os discentes surdos entendem mais através do aspecto visual, levando 
isso em consideração, procurou-se elaborar um material com a menor quantidade de 
textos possível e maior quantidade de elementos visuais, que neste caso foram: imagens, 
diagramas, esquemas e fluxogramas. 

Cada assunto que foi trabalhado em nossos resumos adaptados, foi feito uma 
explicação detalhada do conteúdo abordado, assim descrevemos as variáveis envolvidas 
em cada fórmula e ainda utilizamos imagens nos exemplos resolvidos, para facilitar a 


compreensão e a abstração do conteúdo, conforme mostramos através da figura 1. 
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Cinemática 
E o campo da Física que estuda o movimento com corpos desconsiderando suas 


causas. 
& Ponto material e corpo extenso 


Quando se estuda o movimento de um corpo é necessário saber o comprimento, 
largura e a altura Pois acontecem casos em que as dimensões do corpo são bastante 
pequenas em relação às outras medidas envolvidas, com isso elas podem ser 
desprezadas e, portanto o corpo é considerado um “ponto material” conforme o exemplo 
apresentado na figura 1. 





Fonte: https: “avr coladaweb com fisica mecanica mercia 


Figura 1: Abordagem teórica de um tópico da Física com a associação de imagens 


Fonte: http://coladaweb.com/fisica/mecanica/inercia 


Foram inseridos elementos que representavam as variáveis nas fórmulas de cada 
assunto abordado, também foi feito uma contextualização dos exemplos e exercícios 
resolvidos com imagens, para que não precisasse de muitas informações textuais para 
explicar os assuntos da Física. Cada conteúdo que trabalhamos em nossos resumos 
adaptados, buscamos explicar sempre inserindo variáveis relacionadas com as fórmulas 
da Física, conforme é mostrado na figura 2. 
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Podemos ves a fiumção bosaria da fimção s= f(f) pao MET. 


Fina 2: Faação s= Es) 
= 


2H 





Fembe: Agpositta de sia 9 200 — fincireda 
o O caro parta na preição 5, Po instanta t= D 
Ko seguinte moemento & caro encontrasse na posição final 5 


Ászim partimos da definição da velcridado môdia, que a: 
E 





Dispomos aplicamos a: dafinições acima: 
p= 
Ea 


simplificando a exposssão 
r= TT mp ri=E—sg = Sotmti=s 


I=vlando a posição final na fosmula: 
at mi=s 


> Eo+ mi: 













s = espoço fFinai do corro 
e So SS SSPoço irlciai dor corto 
rr — reslocidoda do COTTo 


ft = tempo do corro 


Figura 2: Contextualização do tema MRU e detalhamento das variáveis envolvidas 


Fonte: Acervo pessoal 


Dessa maneira, procuramos apresentar uma visão mais generalizada do conteúdo, 
só que de maneira simplificada, através da visualização daquilo que foi envolvido como 
variável nos exemplos e exercícios resolvidos, fornecendo as informações necessárias 


para a compreensão do assunto abordado por parte desses discentes surdos. 


Passo 3: Pesquisa e elaboração dos Jogos Educativos Adaptados 


O último passo foi o mais difícil de ser executado, devido à dificuldade de se 
encontrar jogos adaptados sobre temas da Física, encontrou-se vários jogos adaptados 
para surdos, porém todos os jogos que foram encontrados eram voltados para o ensino 
fundamental ou primário, além não encontrarmos muitos sinais em Libras relacionados 
aos termos da Física. No entanto, após realizarmos diversas pesquisas feitas na internet, 
de artigos, monografias, livros digitais, apostilas, bem como o próprio acervo do NAPNE, 


conseguimos um excelente material que nos auxiliou neste trabalho. O material intitulado 
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“Sinalizando a Física” foi a nossa fonte principal dos sinais em Libras, pois ele abrange 
quase todos os termos que são ministrados no 1º ano do Ensino Médio. 

Com o auxílio deste material e da apostila construindo a Física, elaboramos vários 
jogos educativos adaptados para este público como a cruzadinha, o caça-palavras e o jogo 
da memória sobre temas da Física contextualizados com as situações cotidianas destes 


discentes surdos, conforme mostramos nas figuras 3 e 4. 


Crozadinha sobre a Cioematica 


Cinamatica & o estudo do movimento dos * (dica: a resposta tem 
6 letras), 

O 4x representa a variação de + (dica: a resposta tem 6 letras). 
Considare dois observadores À e E em que, o observador À encontra-se dentro 
da um ônibus que se movimenta numa estrada e o observador 6, na beira da 
estrada. Ambos obsarvam uma lâmpada no tato do ônibus, para 0 observador Ea 
lampada esta em movimento e para 0 cbservador À, esta em 7 (dica: 
a reposta tem 7 letras). 

Supondo qua uma formiga estaja atravessando uma rua, nesta caso, ela & 
considerada um ponto matarial ou um corpo extenso” (dica: a resposta tem 13 
letras). 

Uma balsia da 5 metros se dasloca em uma piscina de 12 metros, o referencial 
dala a da ponto matsrial ou um corpo extenso” (dica: a resposta tem 12 letras). 
Qual a orbita que a Tarra se encontra como um ponto material” (dica: a 
resposta tem 11 letras). 


Um aluno esta indo para a aula de carro, em certo momanto o carro faz uma 


parada inesperada, o que acontece com o corpo do garoto” (dica: a resposta 


tem 9 letras), 

Um jogador da selação brasilaira esta cobrando um pênalti a um fã da 
arquibancada da lateral observa a trajetória da bola saindo do chão e chegando a 
trava, qual tipo de trajatória a bola faz! (dica: a resposta tem 9 letras), 





Figura 3: Jogo a Cruzadinha sobre a Cinemática com as dicas das respostas. 


Fonte: Acervo pessoal 
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Figura 4: Jogo sobre a Cinemática com Libras 


Fonte: Acervo pessoal 


Outro jogo educativo adaptado que foi confeccionado com temas da Física foi o 
jogo caça-palavras, que tinha o objetivo de verificar se os discentes tinham assimilado o 
conteúdo que foi trabalhado anteriormente na sala de aula, conforme é mostrado na figura 
5. 
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Caça-palavras sobre as considerações conceituais das Leis de Newton 


As Leis de Newton é a parte da física que aborda o estudo das forças que atuam 
sobre um objeto produzindo ou alterando o seu movimento ou mesmo o deformando-o. As 
Leis de Newton também são conhecidas como os três principios da dinâmica que agem na 
natureza dos corpos, elas são divididas em três leis, mais conhecidas como as Leis de 
Newton, são elas: Princípio da inércia ou primeira lei de Newton; Princípio Fundamental da 
Dinâmica ou segunda lei de Newton; Princípio da ação e reação ou terceira lei de Newton. 

De acordo com o texto sobre as Leis de Newton encontre as palavras em negnto no 


caça-palavras abaixo: 


Lei NE CNEBEHTIMSMO Três 


e 


Inercia 


8 FOEMAMNEIROGI a 
A | 
df IO HVOETTCEABIT 2. 
É | 
Movimento RO E T I E AU I] SA | ta 
Dinâmica 
CERVEIRKAETÕO 
E dã IA EMHROTADOTER| 
im meme | A dês ni nt E 
E* SREAIRLLEONWT|tttivias 


Objeto Newton 


E - 4 Tod. pm [A | e é fil [ = PP foi 

I 1 - a E Ei a! ah E E ai Ra) I E na Pi gi EE fes , E aTái É, h T: 

Wi pr def E E W “q + E =p sê E +] iz f id f Ui] E » ri |? Tr a | E, To 
Terceira Primeira Segunda 


Figura 5: Exemplo do jogo caça-palavras 


Fonte: Acervo pessoal 


Estes jogos educativos adaptados foram propostos após ter sido feita a explanação 
de cada assunto da Física, para que desta forma os discentes surdos não tivessem 
dificuldade e pudessem ter mais autonomia, de acordo com o representado na figura 6. 
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Figura 6: Discentes surdos realizando a atividade do caça-palavras, cujo tema era Leis de 
Newton 


Fonte: Acervo pessoal 


Outro jogo educativo adaptado confeccionado durante o desenvolvimento deste 
projeto, foi justamente o jogo da memória que envolveu dois temas da Física como a 
cinemática e as leis de Newton. Devemos destacar que neste jogo, fio possível fazer uma 
associação de uma imagem que represente o fenômeno da Física com um sinal em Libras 
do termo da Física ou em datilologia, onde é feita a sinalização do alfabeto manual da 


palavra, conforme mostramos na figura 7. 
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Figura 7: Jogo da memória com exemplos de associação de imagens com o sinal de um termo 
da Física em Libras ou em datilologia 


Fonte: Acervo pessoal 


Portanto, os resumos adaptados e os jogos educativos também adaptados que 
foram produzidos, foram disponibilizados no segundo semestre na sala de atendimento 
aos alunos surdos do 1º ano do Ensino Médio no IFAM-CMC, para que estes possam 
servir de material didático de apoio no componente curricular Física e possam ser utilizados 
por outros professores que atendem estes discentes surdos ou outros que adentrarão na 
instituição com o intuito de facilitar a aprendizagem dos mesmos. 

Outra questão a ser considerada é que estes materiais adaptados, também servirão 
de apoio para os professores de Física, que forem ensinar para alunos surdos em suas 
turmas e que muitas vezes não sabem como ensinar para este público específico. Assim, 
este estudante poderá mostrar através dos resumos adaptados e dos jogos, onde está a 
sua dúvida, portanto o aluno poderá pedir explicação do assunto ou na parte onde não está 
entendendo e assim o professor poderá observar como o aluno surdo aprende e poderão 
servir também de modelo caso o professor queira adaptar outro assunto da Física que não 
foi contemplado por este projeto. 
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Passo 4: Acompanhamento das atividades dos discentes surdos nas aulas de 
reforço de Física 

Durante a execução deste projeto foi feito o acompanhamento dos discentes 
surdos na sala de atendimento ao longo de várias aulas de reforço de Física. A sala de 
atendimento é um espaço reservado para realizar as aulas de reforços de várias disciplinas 
como Matemática, Física, Química, Informática, Língua Portuguesa, além da Libras e 
desenvolver atividades extracurriculares com os discentes surdos. 

Foram desenvolvidas várias atividades durante o acompanhamento desses 
discentes surdos, como apoio nas aulas de aulas de reforço sobre temas da Física, auxílio 
na solução de problemas de Física, acompanhamento durante a solução das atividades de 
Física, inclusive durante a aplicação dos jogos educativos adaptados, onde auxiliávamos, 
tiramos as dúvidas dos discentes de acordo com a figura 8. 





Figura 8: Interação dos discentes nas aulas de reforço. 


Fonte: Acervo pessoal 


No decorrer das atividades desenvolvidas por este projeto, percebeu-se a importância 
de se confeccionar materiais adaptados para este público específico, que possui uma forma 


diferente de aprender. 


5 | CONSIDERAÇÕES FINAIS 


O projeto em questão não tem a intensão de fornecer todas as soluções para os 
problemas que rondam a aprendizagem de alunos surdos no Ensino Médio. No entanto, 
iniciativas como esta, aliadas à vontade e ao interesse dos envolvidos, são o que contribui 
para o progresso da educação inclusiva no Brasil e no mundo. Ainda há muito o se que 
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fazer com relação à criação de novos recursos, à capacitação de profissionais e à garantia 
de acesso à educação para as minorias (FERRAZ, 2014). Apesar das dificuldades para 
encontrar os materiais adaptados que servissem de modelo para nos basear, o trabalho 
teve um ótimo resultado, pois a produção dos resumos adaptados e dos jogos educativos, 
auxiliou muito no processo de ensino-aprendizagem desses discentes surdos que 
participaram das atividades do projeto. 

Através dessa iniciativa, percebemos a importância de se confeccionar materiais 
adaptados que possam atender a esse público e que ações como essa devem ser mais 
incentivadas, pois, não se pode deixar esse público a margem do conhecimento, pois 
eles têm os mesmos direitos ao conhecimento que os demais, entretanto, não podemos 
“fechar os olhos” para a situação que está à frente de nossos olhos, pois sabemos que 
este público necessita de uma atenção especial, e que nem sempre é dada, e recursos 
como os que foram confeccionados serão em um futuro não muito distante de fundamental 
importância para auxiliar no processo de ensino-aprendizagem dos alunos surdos, que 
estão adentrando cada vez mais nas escolas regulares ao longo dos anos. 

Com o apoio do NAPNE, que nos concedeu a oportunidade de acompanhar os 
discentes surdos na sala de atendimento, pudemos compreender como estes discentes 
surdos aprendem e desta forma pudemos criar um material que realmente pudesse atender 
o publico em questão, e buscamos produzir com todos os detalhes ao qual deve ser levado 
em consideração, como o fato de que eles entendem mais através do aspecto visual, 
levando isso em consideração procuramos trabalhar com a menor quantidade de textos 
possível e maior quantidade de elementos visuais, que neste caso foram: imagens, tabelas, 
diagramas, esquemas, fluxogramas. Desta forma estes materiais produzidos no decorrer 
deste projeto, como os resumos adaptados e os jogos educativos inclusivos realmente 
contribuíram para a compreensão dos assuntos por parte dos discentes surdos na disciplina 


Física e isso contribuiu para aumentar o rendimento deles na matéria Física. 
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